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1  Introduction 


Figure  1:  Seiko/Epson  pFR  [1] 

1 . 1  Definition  and  Purpose  of  Micro  Air  Vehicles 

Micro  Air  Vehicles,  or  “MAVs”,  are  small  flying  vehicles  that  are  six  inches 
or  less  in  size  and  carry  a  sensing  payload.  An  example  of  a  cutting  edge  MAV  is 
seen  in  Figure  1,  the  Seiko/Epson  “Micro  Flying  Robot”  [1],  MAV  development  has 
been  driven  in  recent  years  primarily  by  military  needs,  with  the  Defense  Advanced 
Research  Project  Agency  (DARPA)  providing  most  of  the  initial  research  impetus.  In 
1996,  DARPA  put  forward  its  original  goals  and  general  definition  for  an  MAV. 
Dimensionally,  the  vehicle  should  be  constrained  to  15  cm  with  a  weight  of  about  100 
grams.  The  endurance  goal  was  encouraged  to  be  one  hour  and  the  payload  should 
include  a  camera  or  other  sensing  device.  Today,  however,  some  vehicles  with  a 
weight  of  300  grams  or  more  are  still  considered  MAVs.  No  published  MAV  design 
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with  a  size  of  six  inches  or  smaller  has  apparently  reached  the  goal  of  one  hour 
endurance. 

For  the  military,  an  MAV  can  provide  stealthy  surveillance  into  a  complex 
and  possibly  dangerous  situation  without  much  risk.  Common  environments  for 
usage  include  urban  areas  and  inside  buildings  or  in  complex  terrain  such  as  hills, 
mountains  or  inside  caves.  For  civilian  applications,  an  MAV  can  examine  an 
environment  that  is  harmful  due  to  structural,  chemical,  electrical,  or  other  hazardous 
concerns.  An  MAV  could  carry  a  video  camera  or  an  array  of  other  sensors,  such  as 
those  to  sense  heat,  chemicals,  or  radioactivity.  The  goal  for  portability  of  MAVs  is 
for  one  person  to  be  able  to  carry  the  entire  system.  All  of  the  components,  including 
ground  based  instruments,  should  be  able  to  be  quickly  assembled  and  disassembled 
and  packed  into  a  backpack.  It  is  expected  that  an  MAV  should  only  require  one  or 
two  people  to  transport,  assemble,  and  operate. 

1.2  State  of  the  Art:  Current  MAV  configurations 

While  there  is  no  perfect  MAV  for  all  situations,  most  of  the  small  flying 
vehicles  being  developed  can  be  divided  into  three  categories.  Just  like  full  size 
aircraft,  the  two  most  successful  configurations  are  fixed  wing  and  rotary  wing.  The 
third  configuration,  which  imbibes  nature,  is  a  flapping  wing  configuration.  While 
functional  fixed  wing  and  rotary  wing  MAVs  are  numerous  and  successful,  there  are 
only  a  few  experimental  flapping  wing  MAVs  that  have  not  been  used  in  the  field. 
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1.2.1  Fixed  Wing 


Fixed  wing  MAVs  have  so  far  been  the  most  successful  at  achieving  the 
longest  endurance  and  greatest  speed  and  range.  They  are  relatively  simple,  fast,  and 
efficient  compared  to  other  categories  of  MAVs.  The  BAE/Lockheed  Martin 
Microstar,  seen  in  Figure  2,  is  an  example  of  a  typical  fixed  wing  MAV.  It 
incorporates  a  part  of  its  support  systems  into  its  structure,  such  as  the  patented  circuit 
board  that  is  built  as  part  of  the  wing  [2],  Due  to  many  innovations  to  save  weight,  it 
has  an  endurance  of  about  25  minutes  and  can  fly  at  speeds  of  30  kts  while  carrying  a 
video  camera  payload. 


Figure  2:  BAE/Lockheed  Martin  Microstar  MAV 


The  Aerovironment  Black  Widow,  seen  in  Figure  3,  also  incorporates  its 
systems  and  payload  within  the  wing,  providing  a  large  wing  area  to  lift  its  mass  of 
80  grams.  It  has  an  endurance  of  up  to  30  minutes  as  a  result  of  optimization  of  its 
wing  planform  and  cross  section,  electrical  systems,  and  direct  drive  motor  and 
propeller.  During  testing,  a  larger  propeller  was  found  more  aerodynamically 
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efficient  when  paired  with  a  combination  motor  and  gearbox.  The  gearbox 
efficiency,  however,  is  only  81%  because  of  its  small  size,  which  negated  the  increase 
in  efficiency  of  a  larger  diameter  propeller  [3].  These  fixed  wing  craft  have 
constraint  on  their  minimum  airspeed,  which  is  not  conducive  to  flying  indoors  or  in 
confined  spaces.  Because  of  their  high  minimum  airspeed,  they  are  not  capable  of 
“hover  and  stare”,  “perch  and  stare”,  or  tagging  and  targeting  operations. 


Figure  3:  Aerovironment  “Black  Widow” 


1.2.2  Rotary  Wing 

For  very  low  speed  and  hover  capability,  rotary  wing  MAVs  are  most 
frequently  utilized.  Just  like  full  size  helicopters,  these  vehicles  use  either  a  single 
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main  rotor  with  a  mechanism  for  anti-torque  or  a  coaxial  configuration.  In  the  MAV 
realm,  coaxial  designs  appear  attractive  because  of  their  inherent  compactness. 
Examples  of  coaxial  MAVs  are  the  Seiko  Micro  Flying  Robot  and  the  Micor, 
developed  at  the  University  of  Maryland  as  seen  in  Figure  4  [4],  The  Micor’ s  rotors 
are  counter-rotating  with  an  equal  and  opposite  torque  applied  to  each  rotor  to  lift  the 
vehicle’s  150  gram  mass.  The  vehicle  is  controlled  only  by  the  teetering  upper  rotor. 
Through  innovations  in  low  Reynolds  number  blade  design  and  optimized  drivetrain, 
the  vehicle  is  able  to  hover  for  up  to  20  minutes.  The  disadvantages  of  this  design 
include  aerodynamic  interference  between  the  two  rotors  and  mechanical  complexity 
and  added  weight  of  a  dual  rotor  system. 


Unlike  a  coaxial  design,  single  main  rotor  designs  require  a  separate  device 
for  anti-torque.  Some  concepts  are  basically  scaled  down  conventional  helicopters, 
with  a  large  main  rotor  and  a  small  tail  rotor,  such  as  the  radio  controlled  “Micron” 
helicopter  seen  in  Figure  5  [5].  The  57  gram  Precision  Heli  Micron  V2  has  an  8” 
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fixed  pitch  rotor  that  controls  vehicle  thrust  by  varying  the  rotor  RPM  and  an 
endurance  of  about  10  minutes.  The  benefits  of  this  conventional  design  are 
simplicity  and  familiarity  in  construction  and  control.  The  tail  rotor  and  boom, 
however,  add  to  the  vehicle  dimensions  and  are  a  detriment  to  the  goal  of 
compactness.  Scaling  down  a  conventional  helicopter  design  degrades  the  efficiency 
of  the  vehicle.  While  the  main  rotor  of  a  conventional  helicopter  operates  at  a  figure 
of  merit  as  high  as  0.85,  optimized  MAV  rotors  currently  operate  with  a  figure  of 
merit  of  about  0.6  [6]. 


Figure  5:  Precision  Heli  Micron  V2 


In  order  to  meet  the  goal  of  a  compact  and  simple  design  without  the  need  for 
two  rotors,  some  MAVs  utilize  a  single  rotor  with  vanes  in  the  downwash  to 
counteract  the  torque  required  by  the  main  rotor.  As  the  main  rotor  torque  and  thrust 
increase,  the  anti-torque  from  the  vanes  also  increases  due  to  the  increased 
downwash.  Three  examples  of  vanes  used  as  a  stable  method  of  control  and  torque 
counteraction  are  the  Honeywell  iSTAR  [7]  and  the  Giant  and  the  TiShrov  at  the 
University  of  Maryland  [8,  9,  10].  The  Giant  is  seen  in  Figure  6  and  the  TiShrov  in 
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Figure  7.  The  Giant  and  the  TiShrov  are  similar  vehicles,  but  the  TiShrov  utilizes  a 
shrouded  rotor.  The  shroud  has  a  curved  inlet  that  increases  aerodynamic  efficiency 
by  reducing  tip  losses  and  providing  additional  thrust  due  to  the  accelerated  airflow 
over  the  inlet.  The  thrust  benefits  of  the  shroud  have  not  yet  exceeded  the  shroud’s 
weight,  so  the  unshrouded  Giant  has  a  greater  endurance  of  about  8  minutes. 


Figure  6:  University  of  Maryland  “Giant” 


The  Giant  incorporates  a  fixed  pitch  rotor  with  a  diameter  of  135  cm  and  a 
Hiller  stabilizer  bar  for  control.  Two  of  its  vanes  have  articulated  flaps  to  trim  the 
vehicle  with  anti-torque  and  provide  yaw  control.  The  vanes  on  the  Giant  and 
TiShrov  have  been  tested  to  consume  approximately  10%  of  the  vehicle’s  power, 
which  is  comparable  to  the  tail  rotor  on  a  full  scale  helicopter.  The  disadvantage  of 
the  current  design  of  the  Giant’s  vanes  and  the  TiShrov’s  shroud  is  that  they  cause 
more  drag  than  a  conventional  rotor. 
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Figure  7:  University  of  Maryland  “TiShrov” 


The  Honeywell/Allied  Aerospace  iSTAR  is  seen  in  Figure  8.  Like  the 
TiShrov,  the  iSTAR  uses  a  shroud  and  vanes  in  the  downwash  of  its  single  main  rotor 
to  counteract  the  rotor  torque  and  also  for  pitch  and  yaw  control.  At  most,  the  iSTAR 
has  been  tested  with  a  29”  rotor  and  duct  to  increase  payload  and  endurance.  While 
most  rotary  wing  MAVs  can  only  be  operated  at  slow  speeds  and  minimal  wind,  the 
iSTAR  design  has  the  capability  to  transition  to  high  speed  forward  flight.  The 
vehicle  can  tilt  almost  90°  and  gain  horizontal  airspeed  so  that  its  rotor  acts  like  a 
ducted  propeller  and  the  duct  acts  as  a  circular  wing  to  provide  lift.  In  order  to 
achieve  the  amount  of  control  authority  necessary  to  pitch  the  vehicle  for  horizontal 
axial  flight,  the  vehicle  has  a  high  vertical  center  of  gravity.  The  vehicles  engine  can 
be  seen  in  Figure  8  mounted  above  the  main  rotor  to  raise  the  center  of  gravity.  The 
high  center  of  gravity  leads  to  a  difficulty  in  pitch  control,  especially  at  intermediate 
forward  flight  speeds  [11],  The  vehicle  also  has  a  high  profile  drag  from  its  shroud 
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and  a  loud  noise  signature  from  its  piston  engine. 


Figure  8:  iSTAR  at  rest  and  transitioning  to  forward  flight 
1.2.3  Flapping  Wing 

In  addition  to  the  fixed  wing  and  rotary  wing  designs  for  MAVs,  there  are  also 
some  unconventional  designs  that  utilize  flapping  wings  to  generate  lift  and 
propulsion.  Most  of  these  concepts  are  based  upon  nature  and  simulate  the  motions 
of  the  wings  of  flying  creatures.  The  two  basic  types  of  flapping  motion  found  in 
nature  are  avian  based  and  insect  based.  In  pure  avian  based  flapping,  the  wings 
stroke  in  a  nearly  vertical  plane  with  respect  to  the  body  frame.  The  wing  of  a  bird  is 
intricately  woven  with  muscles  for  controlling  the  pitch,  camber,  and  sweep  of  the 
wing  as  well  as  controlling  the  orientation  of  individual  feathers  depending  on  the 
flight  condition.  In  this  regard,  the  performance  of  avian  flight  is  extremely  reliant  on 
active  morphing  for  maximum  flight  efficiency.  For  most  birds  and  avian  based 
flight  vehicles,  hovering  is  not  possible.  They  usually  depend  on  a  forward  airspeed 
to  generate  lift  and  utilize  flapping  to  produce  a  forward  thrust.  The  Aerovironment 
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Microbat,  seen  in  Figure  9,  is  only  one  example  of  an  MAV  whose  flapping  motion  is 
based  upon  avian  flight  [12].  The  wings  of  the  Microbat  do  not  actively  change 
shape,  however.  They  change  shape  passively  due  to  inertial  and  aerodynamic  forces. 
The  wings  flap  simply  in  a  vertical  plane  to  produce  thrust  and  lift.  Like  most  birds 
of  its  size,  the  Microbat  requires  a  horizontal  velocity  to  fly  and  it  cannot  hover. 


Figure  9:  Aerovironment  “Microbat” 


In  hover  capable  insect  based  flapping,  the  stroke  plane  is  nearly  horizontal, 
with  the  wing  changing  direction  with  respect  to  the  flowfield  from  the  upstroke  to 
the  downstroke,  as  seen  in  Figure  10  [13].  The  wing  pronates  and  supinates,  or 
inverts,  each  time  it  changes  directions.  The  wings  create  their  own  freestream 
velocities  and  their  net  force  averages  to  a  vertical  force  during  hover  when  the  stroke 
plane  is  horizontal.  In  forward  flight,  the  stroke  plane  tilts  forward,  which  tilts  the  net 
thrust  vector  forward.  Like  avian  based  flapping,  the  wings  change  shape  depending 
on  flow  conditions.  Unlike  avian  based  flapping,  however,  changes  in  the  wing  shape 
occur  purely  passively  due  to  aeroelastic  coupling  of  inertial  and  aerodynamic  forces. 
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Net  Force 


Figure  10:  Insect  based  flapping 


The  Mentor  [14],  seen  in  Figure  11,  fits  the  description  of  insect  based  flight. 
With  a  wingspan  of  12”  and  an  endurance  of  about  10  minutes,  it  was  the  first  small 
flapping  wing  vehicle  capable  of  hovering  and  transitioning  into  forward  flight.  The 
four  wings  flap  in  a  horizontal  stroke  plane  and  use  the  “clap-fling”  phenomenon 
presented  by  Weis-fogh  [15].  Clap-fling  occurs  on  the  Mentor  when  the  wings  are 
forced  together  and  their  leading  edges  make  contact  first.  The  flexible  wing 
membranes  squeeze  the  air  in  between  them,  forming  an  area  of  high  pressure  which 
generates  lift.  The  flexible  wings  allow  the  higher  pressure  air  to  escape  downwards 
as  they  clap  together.  When  the  wings  are  pulled  apart,  the  flexible  membranes  again 
lag  the  motion  of  the  stiffer  leading  edge,  forming  an  area  of  low  pressure  on  their 
upper  surface.  This  low  pressure  region  creates  lift  and  also  pulls  air  downwards, 
generating  momentum  flow.  The  negative  effects  of  this  configuration  include  short 
endurance  and  a  loud  noise  signature  and  material  fatigue  due  to  the  high  speed 
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clapping  of  the  wings.  None  of  the  flapping  wing  MAVs  have  been  successful  at 
achieving  the  endurance  and  range  of  current  fixed  wing  or  rotary  wing  designs. 


Figure  11:  “Mentor”  flapping  wing  vehicle 


1.3  MAV  Challenges 

During  the  past  decade,  electronic  components  have  become  more 
miniaturized  and  lightweight.  Only  recently  have  they  become  feasible  for  MAV  use. 
Even  with  these  recent  advancements,  no  published  MAV  with  a  size  of  six  inches  or 
less  has  reached  the  endurance  goal  of  one  hour  sustained  flight.  One  cause  of  the 
endurance  shortfall  is  the  weight  of  the  payload  and  electronics  for  navigation  and 
control.  As  electronics  and  MEMS  technologies  are  further  miniaturized,  MAV 
payloads  will  be  manufactured  lighter,  allowing  higher  performance,  reduced  size,  or 
increased  endurance.  Another  cause  for  short  endurance  is  the  size  and  efficiency  of 
the  power  source  and  motor.  While  large  electric  motors  have  efficiencies  greater 
than  90%,  MAV  sized  motors  have  efficiencies  on  the  order  of  60%  or  less,  which 
limit  performance  [3].  As  the  motor  efficiency  and  the  battery  specific  energy  and 
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specific  power  increase,  there  will  be  an  improvement  in  endurance.  Also  to  improve 
endurance,  there  is  an  opportunity  to  optimize  the  aerodynamic  design  of  MAVs  for 
operation  in  a  low  Reynolds  number  regime.  Many  of  the  aerodynamic  theories  and 
assumptions  applied  to  full  scale  aircraft,  such  as  turbulent  boundary  layers  and 
steady,  attached  flows,  may  no  longer  be  valid  as  the  scale  of  flying  vehicles  becomes 
small. 

1.3.1  Low  Reynolds  Number  Aerodynamics 

A  defining  characteristic  of  low  Reynolds  number  aerodynamics  is  laminar 
flow.  As  the  thrusting  and  lifting  surfaces  of  a  flying  vehicle  become  smaller,  they 
become  inherently  less  efficient.  Compared  to  the  pressure  forces  and  inertial  forces 
on  the  vehicle,  the  relative  fluidic  viscous  forces  become  greater  as  the  vehicle  size 
decreases.  The  ratio  of  vehicle  inertial  forces  to  viscous  forces  is  commonly 
expressed  in  terms  of  Reynolds  number.  Reynolds  number  is  calculated  in  terms  of 
fluid  density  (p ),  free-stream  velocity  (v),  a  reference  length  of  the  object  (L),  which 
is  usually  chord  length  of  a  wing,  and  a  dynamic  viscosity  constant  of  the  fluid  (/./): 

„  pvL  Inertial  force 

Re  =  - —  = - 

//  Viscous  force 

As  seen  in  Figure  12,  the  Reynolds  Numbers  for  MAVs  can  be  on  the  order  of 
10,000  or  even  lower  [4],  Many  previous  experiments  have  shown  that  airfoils 
optimized  for  full  scale  aircraft  and  helicopters  perform  poorly  at  MAV  Reynolds 
numbers.  The  detrimental  effects  of  low  Reynolds  numbers  include  reduced 
maximum  lift  coefficients  and  increased  drag  coefficients  due  to  laminar  flows  and 
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flow  separation.  These  effects  have  been  measured  in  wind  tunnel  experiments  by 
Laitone  [16],  Schmitz  [17],  Sunada  et  al  [18],  Selig  et  al  [19],  Mueller  [20,  21]  and 
many  others. 


Re 

Figure  12:  Vehicle  Mass  vs.  Reynolds  Number 


The  results  show  when  moving  from  Reynolds  numbers  on  the  order  of  106 
down  to  Reynolds  numbers  of  20,000  and  below,  the  maximum  lift  coefficients  of 
conventional  airfoils  can  be  reduced  by  a  factor  of  two.  Not  only  does  the  lift 
performance  degrade,  but  the  drag  coefficients  are  many  times  higher.  This  causes 
maximum  lift  to  drag  ratios  of  conventional  airfoils  at  low  Reynolds  numbers  to  be  an 
order  of  magnitude  less  than  their  values  at  higher  Reynolds  numbers.  Increased 
profile  drag  and  induced  drag  are  both  responsible  for  increased  drag  at  low  Reynolds 
numbers.  Due  to  these  effects,  an  MAV  rotor  requires  a  profile  power  that  is  50% 
higher  than  a  full  scale  helicopter  rotor  and  suffers  from  lift  to  drag  ratios  ranging 
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from  about  2  to  8  [22], 

Figure  13  shows  well  known  airfoil  data  in  the  form  of  Cl  vs.  a  and  Cd  vs.  a 
for  a  NACA  4412  airfoil  that  produces  very  high  lift  to  drag  ratios  at  full  scale  [23, 
24],  It  is  obvious  that,  at  Reynolds  numbers  of  20,000,  the  performance  of  this  airfoil 
is  drastically  different  from  its  behavior  at  higher  Reynolds  numbers.  The  lift  to  drag 
ratio  is  about  100  at  Reynolds  numbers  of  3  x  106,  but  falls  to  about  5  at  a  Reynolds 
number  of  20,000.  This  is  a  typical  MAV  lift  to  drag  ratio  and  it  causes  one  of  the 
biggest  challenges  of  applying  conventional  lift  and  propulsion  methods  to  MAVs. 


Figure  13:  Cl  vs.  a  and  Cl  vs.  Cd  at  different  Reynolds  numbers  [23,  24] 


Similar  Reynolds  number  effects  were  found  by  Schmitz  [17],  who  performed 
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many  tests  of  common  airfoils  at  Reynolds  numbers  ranging  from  20,000  to  160,000. 
It  was  found  that,  decreasing  through  this  Reynolds  number  range,  maximum  lift 
coefficients  can  fall  by  a  factor  of  more  than  two  and  minimum  drag  coefficients  can 
double.  A  “critical  Reynolds  number”  was  identified,  which  is  on  the  order  of  60,000 
to  100,000,  depending  on  the  airfoil.  Above  this  critical  Reynolds  number,  the  flow 
over  the  airfoil  was  found  turbulent  and  the  skin  friction  drag  was  relatively  low. 
There  was  laminar  flow  below  the  critical  Reynolds  number,  which  reduced 
maximum  lift  coefficients  and  increased  skin  friction  drag  dramatically.  Formulas 
were  presented  to  predict  how  drag  coefficients  increase  and  lift  coefficients  decrease 
as  Reynolds  number  decreases. 

Researchers  have  found  that  the  lift  curves  of  flat  plates  are  relatively 
unaffected  by  low  Reynolds  numbers.  While  conventional  airfoils  outperform  thin 
flat  plates  at  high  Reynolds  numbers,  the  plates  do  not  suffer  the  two  or  threefold 
reduction  in  maximum  lift  coefficients  below  their  critical  Reynolds  number.  Figure 
14  shows  the  lift  vs.  drag  performance  of  a  conventional  airfoil,  a  flat  plate,  and  a 
cambered  plate  above  and  below  the  critical  Reynolds  number  [25],  At  low  Reynolds 
numbers,  conventional  airfoils  will  have  laminar  flow  along  their  surface,  which 
separates  readily  and  stalls  the  airfoil  at  low  angles  of  attack.  The  abrupt  and 
sometimes  sharpened  leading  edge  of  a  flat  or  curved  arc  airfoil  acts  as  a  boundary 
layer  trip,  causing  the  flow  to  become  turbulent  once  it  passes  the  leading  edge  of  the 
airfoil.  Even  with  laminar  flow,  the  lifting  performance  of  flat  plate  airfoils  is  not 
significantly  affected  by  passing  below  the  critical  Reynolds  number.  This  is  the 
reason  that  many  current  MAVs  utilize  curved  plates  or  membrane  wings. 
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R  =  40,000 


R  =  120,000 


Figure  14:  Cl  vs.  Cd  for  three  airfoils 
at  Reynolds  numbers  of  40,000  and  120,000  [25] 

1 .4  Previous  Research 

Whether  applying  low  Reynolds  number  airfoil  designs  to  small  helicopter 
rotors  or  studying  phenomena  from  nature  for  lift  and  propulsion,  the  previous 
research  presented  had  a  common  goal  of  optimizing  and  understanding  small  scale 
flight.  The  nature  based  research  fits  into  two  categories.  The  first  category  is 
“biomimetic”,  which  means  imitating  the  motions  and  planforms  found  in  nature. 
The  second  category  is  “bio-inspired”  and  it  usually  utilizes  oscillating  wings  in  a 
slightly  different  manner  than  in  nature. 
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1 .4. 1  Conventional  Micro  Rotors 


It  is  now  well  established  that  optimal  designs  for  conventional  airfoils  for  full 
scale  aircraft  are  not  optimal  designs  for  MAVs.  Recent  research  into  micro  rotor 
performance  has  verified  past  findings  of  low  Reynolds  numbers  airfoil  designs, 
which  show  that  the  performance  of  thin  cambered  plates  exceeds  that  of 
conventional  airfoils  at  low  Reynolds  numbers  [16,  17,  19].  With  a  decrease  in 
Reynolds  numbers,  the  lift  coefficients  of  thin  plates  are  relatively  unchanged,  while 
the  lift  capability  of  conventional  airfoils  is  reduced  significantly. 

At  the  University  of  Maryland,  there  has  been  much  research  accomplished  on 
aerodynamic  performance  improvements  and  optimization  of  airfoils  and  blade 
planforms  for  MAV  rotors.  Bohorquez  and  Pines  [4]  studied  effects  of  rotor  solidity, 
airfoil  selection,  blade  twist  and  planform  shape  on  MAV  rotor  efficiency  while  using 
simple  and  consistent  manufacturing  methods.  The  study  began  with  an  8.8” 
diameter  two-bladed  rotor  from  a  flying  MAV  spun  at  speeds  between  1333  and  3500 
RPM,  which  resulted  in  tip  Reynolds  numbers  on  the  order  of  20,000  -  40,000.  The 
first  blades  tested  were  rectangular  flat  plates  and  resulted  in  a  figure  of  merit  of  0.35. 
Through  experimental  improvements  of  the  blade  design,  the  final  micro  rotor  had  a 
figure  of  merit  up  to  0.65.  The  experimentally  optimized  blades  incorporated  a  linear 
asymmetric  taper  ratio  of  2:1  starting  at  60%  of  the  rotor  radius  and  gave  the  rotor  a 
solidity  of  0.12.  The  final  airfoils  used  were  2.5%  thick  plates  cambered  by  9%, 
which  is  a  drastic  change  from  the  10-12%  thick  airfoils  commonly  used  on  full  scale 
helicopters. 
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Hein  and  Chopra  [26]  also  focused  on  improving  micro  rotor  performance  at 
low  Reynolds  numbers.  They  studied  different  blade  planforms,  twist  distributions, 
and  airfoil  shapes  on  a  6”  diameter  MAV  rotor  with  tip  Reynolds  numbers  between 
25,000  and  40,000.  The  rotor  performance  was  improved  by  increasing  the  solidity 
above  that  of  full  scale  rotors  to  a  thrust  weighted  solidity  of  0.16.  It  was  found  that 
rotor  efficiency  was  only  marginally  improved  with  moderate  amounts  of  blade  twist 
or  taper.  The  focus  shifted  to  improving  the  airfoil  shapes  used  in  rotary  wing  MAV 
blades.  A  2.75%  thin  rectangular  plate  cambered  7%  with  a  sharpened  leading  edge 
provided  far  superior  results  to  conventional  airfoils  at  this  scale.  Using  this  blade, 
the  test  rotor  achieved  a  figure  of  merit  of  0.55,  which,  while  outstanding  for  low 
Reynolds  numbers,  is  not  close  to  matching  the  efficiency  of  full  scale  rotors.  Flow 
visualization  showed  large  tip  vortices,  a  very  non  uniform  inflow,  including  a 
blocked  flow  near  the  blade  roots,  and  higher  induced  velocities  near  the  blade  tips 
compared  to  full  scale  rotors. 

1 .4.2  Biomimetic  Movement 

When  conventional  steady  aerodynamic  analysis  is  applied  to  the  wings  of  a 
bumblebee,  the  results  show  that  the  bumblebee  wings  cannot  lift  its  weight.  In  order 
for  the  flight  of  many  animals  to  be  possible,  they  achieve  greater  lift  coefficients 
than  can  be  explained  by  conventional  steady  aerodynamics.  Some  researchers  are 
concentrating  on  the  aerodynamics  and  exact  wing  motions  of  animal  flight.  Through 
examinations  of  live  insects  in  flight,  Ellington  [27,  28]  and  Dickinson  [29]  have 
carefully  studied  the  kinematics  of  the  creatures  and  observed  the  different 
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aerodynamic  techniques  used  in  nature  to  achieve  high  vertical  lift  despite  low 
Reynolds  numbers. 

The  simplest  and  most  common  phenomena  used  by  all  hovering  creatures  are 
delayed  stall,  rotational  lift,  and  wake  capture.  A  diagram  from  Dickinson  [30]  in 
Figure  15  shows  how  these  three  aerodynamic  mechanisms  are  used  together  in  insect 
flight.  The  black  line  segment  represents  a  wing  section.  The  red  arrows  show  the 
wing  motion  and  the  black  curved  arrows  show  the  motion  of  the  air  around  the  wing. 
The  blue  arrows  show  the  resultant  force  on  the  wing  section,  which  includes  lift  and 
drag.  This  figure  shows  the  motion  of  a  simple  flat  wing  section  moving  from  right 
to  left  in  steps  1-3  and  then  reversing  direction  and  moving  from  left  to  right  in  steps 
4  and  5.  In  step  1,  when  the  wing  first  starts  translating  at  a  high  angle  of  attack,  the 
wing  does  not  stall.  Instead,  the  pressure  above  the  wing  remains  low  due  to  the 
leading  edge  vortex.  Lower  pressure  at  the  tips  causes  spanwise  flow  that  gradually 
sheds  the  vorticity  of  the  leading  edge  vortex.  Before  the  vortex  can  shed,  the  wing 
translation  slows  and  simultaneously  rotates  in  steps  2  and  3,  as  if  to  a  higher  angle  of 
attack.  Adding  circulation  to  the  flow  generates  lift,  and  the  wing  rotation  in  steps  2 
and  3  generates  more  circulation  than  translation  alone.  Finally,  the  wing  reverses 
direction  and  the  viscous  wake  from  the  previous  stroke  interacts  with  the  wing, 
thereby  increasing  its  lift.  Utilizing  such  mechanisms  in  flight,  insects  are  able  to 
create  much  higher  average  lift  coefficients  than  steady  aerodynamics  might  imply. 
A  fruit  fly,  for  example,  produces  an  average  lift  coefficient  of  about  1.4  over  the 
complete  cycle  of  its  stroke  [30],  Even  the  maximum  lift  coefficients  of  optimized 
airfoils  cannot  match  the  average  lift  coefficients  that  insects  produce  at  this  triple 
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digit  Reynolds  number. 


Delayed  Stall 

Leading-edge  Net  aerodynamic 


Rotational  Lift 


Figure  15:  Delayed  Stall,  Rotational  Lift,  Wake  Capture 


Once  the  kinematic  motions  of  flying  creatures  were  understood,  researchers 
began  simulating  the  motions  with  mechanical  apparatuses  in  the  lab.  Some  hoped  to 
build  a  flying  vehicle,  but  the  most  useful  research  has  been  through  flapping  models 
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used  to  better  understand  the  aerodynamics.  In  these  mechanical  models,  forces  are 
measured  directly  and  the  flow  can  be  carefully  scrutinized.  The  leading  edge  vortex 
that  delays  stall  and  allows  higher  angles  of  attack  has  been  verified  and  shown  to  be 
a  stable  phenomenon.  It  remains  attached  even  under  steady  translation,  because  it 
spirals  out  toward  the  lower  pressure  region  near  the  wingtip,  where  it  combines  with 
the  wingtip  vortex  [31,  32],  The  leading  edge  vortex  remains  attached  due  to 
spanwise  flow.  It  has  been  suggested  through  theory  and  experiment  that  the 
spanwise  flow  and  stable  leading  edge  vortex  occurs  only  over  a  range  of  Reynolds 
numbers  on  the  order  of  10  -10  [31,  33].  MAVs  normally  operate  within  this 
range. 

1.4.3  Bio-Inspired  Aerodynamics 

Some  of  the  same  researchers  that  are  simulating  insect  wing  motion  are  also 
using  simpler  methods  and  models  to  study  the  aerodynamic  phenomena  associated 
with  flapping  motion.  Testing  fixed  wings  of  hover  capable  flying  creatures  in  a 
wind  tunnel  does  not  provide  accurate  results,  because  it  subjects  the  entire  wing  to  a 
uniform  freestream  velocity  [34],  Each  wing  section  of  a  hovering  animal  is  not 
subjected  to  the  same  freestream  velocity  because  the  wing  is  flapped  while  it  is  fixed 
at  the  root.  The  velocities  at  the  tip  will  therefore  be  higher,  just  like  a  helicopter 
blade.  The  steady  aerodynamics  that  occur  during  mid-stroke  can  be  observed  by 
rotating  the  wings  on  a  hover  stand,  visualizing  the  flow,  and  measuring  the  forces. 

Many  rotating  animal  wing  experiments  have  been  performed  by  Usherwood 
and  Ellington  [34,  35]  and  Altshuler,  Dudley  and  Ellington  [36]  at  low  Reynolds 
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numbers  (<  10,000)  with  actual  wings  and  scale  models  of  wings  from  different 
insects  and  birds.  The  researchers  used  a  rotating  hover  stand  because  the  spanwise 
flow  that  stabilizes  the  leading  edge  vortex  is  weak  in  a  conventional  wind  tunnel  test. 
The  rotating  wing  experiments  have  shown  a  stable  leading  edge  vortex  and  spanwise 
flow  during  steady  rotation.  Due  to  the  stable  leading  edge  vortex,  angles  of  attack 
reached  upwards  of  40°  without  excessive  flow  separation  and  stall.  In  fact,  the  lift 
curve  slopes  during  these  experiments  remained  linear  up  to  and  above  30°  angle  of 
attack  and  maximum  lift  coefficients  were  about  1.3  for  most  of  the  wings  tested. 
Actual  hummingbird  wings  showed  the  highest  maximum  lift  to  drag  ratio  of  10, 
while  the  rest  of  the  wings  tested  produced  lift  to  drag  ratios  of  about  5.  Most  insects 
in  hover  achieve  average  lift  coefficients  higher  than  the  lift  coefficients  that  their 
wings  can  produce  in  steady  motion,  whether  rotation  or  translation.  Therefore,  other 
unsteady  mechanisms,  such  as  delayed  stall  and  rotational  lift,  are  clearly  at  work. 

There  have  literally  been  dozens  of  studies  and  experiments  completed  on  the 
subject  of  dynamic  and  delayed  stall  of  translating  and  rotating  wings  and  airfoils. 
Since  most  of  these  studies  were  conducted  to  understand  full  size  helicopters  and 
aircraft,  relatively  few  have  been  done  that  focus  on  low  Reynolds  number  delayed 
stall.  One  low  Reynolds  number  study  by  Okamoto  and  Azuma  focuses  on  heaving 
and  feathering  oscillations  of  wings  in  a  wind  tunnel  at  a  Reynolds  number  of  8000 
[37],  Instantaneous  and  average  forces  were  measured  while  the  wings  were 
operating  both  steadily  and  oscillating  at  a  reduced  frequency  of  0.32  in  heaving, 
feathering,  or  a  combination  of  heaving  and  feathering.  Flat  plates,  cambered  plates, 
and  corrugated  plates  were  tested,  all  of  which  are  known  to  perform  well  at  low 
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Reynolds  numbers. 

The  flat  plates  showed  evidence  of  stall  under  steady  conditions  at  8°  angle  of 
attack  and  a  lift  coefficient  of  0.7.  While  undergoing  heaving  motion  at  a  reduced 
frequency  k  of  0.32,  the  maximum  lift  coefficient  increased  to  1.1  at  an  average  angle 
of  attack  of  26°.  That  is  an  increase  in  average  maximum  lift  coefficient  of  over  50%. 
A  curved  plate  with  9%  camber,  which  generally  achieves  one  of  the  highest  lift 
coefficients  in  this  flow  regime,  was  also  tested.  In  steady  state,  it  achieved  a 
maximum  lift  coefficient  of  1.2  at  an  angle  of  attack  of  12°.  During  periodic  heaving 
motion,  the  maximum  average  lift  coefficient  increased  about  20%  to  1.4.  With  an 
optimal  combination  of  heaving  and  feathering  at  a  phase  difference  of  90°,  these 
wings  produced  positive  thrust  coefficients. 

The  idea  of  plunging  wings  creating  a  propulsive  force  is  not  new.  As  early  as 
1909  and  1912,  Knoller  and  Betz,  respectively,  independently  theorized  that  a  wing 
can  produce  thrust  when  it  plunges  periodically  with  respect  to  a  freestream  velocity 
[38,  39],  They  were  investigating  the  method  for  avian  propulsion.  In  1922,  the 
Knoller-Betz  effect  was  experimentally  verified  by  Katzmeyer  with  wind  tunnel 
measurements  [40],  The  idea  behind  the  effect  is  quasi-steady  in  nature  and  is  seen  in 
the  simplified  diagram  in  Figure  16. 
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/  ^  T 

Figure  16:  Forces  and  velocities  on  a  plunging  airfoil 


The  plunging  airfoil’s  vertical  position  is  defined  by  z(t).  The  vertical 
velocity  of  air  induced  on  the  airfoil  is  in  the  opposite  direction  of  the  airfoil’s 
vertical  velocity  and  is  shown  as  -z(t).  When  the  vertical  air  velocity  -z(t)  combines 
with  the  freestream  velocity  Vx,  the  resultant  freestream  velocity  is  Vefj,  which  causes 
an  effective  angle  of  attack  aeff.  The  lift  L  acts  perpendicular  to  Veff,  which  means  the 
lift  tilts  in  a  forward  direction.  The  lift  force  here  is  broken  down  into  two  forces  on 
the  airfoil:  the  normal  force  N  and  the  thrust  force  T.  If  the  lift  is  great  enough,  the 
plunging  airfoil  can  produce  a  useful  thrust.  The  drag  in  this  figure  is  assumed  small 
with  respect  to  the  lift. 

Recently,  the  Naval  Research  Lab  has  examinined  plunging  and  pitching 
wings  as  a  primary  propulsion  method  for  MAVs.  Their  vehicles,  such  as  the 
“delphinopter”  in  Figure  17,  have  been  some  of  the  most  successful  flying  vehicles  of 
any  scale  that  utilize  this  propulsion  method.  The  delphinopter  produces  thrust  with 
its  plunging  and  pitching  rear  wing  and  produces  lift  with  its  fixed  front  wing  [41]. 
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1 .4.4  Flapping  Rotor  Concept 

If  flapping  wings  can  be  the  propulsion  mechanism  for  birds  and  aircraft,  then 
flapping  blades  should  be  able  to  power  the  rotation  of  a  helicopter  rotor  as  well. 
This  is  the  expectation  of  a  research  group  at  Delft  University  of  Technology  [42], 
They  have  developed  a  flapping  rotor  system  that  flaps  its  blades  once  per  revolution 
at  an  amplitude  of  12°.  The  amplitude  is  enough  for  the  lift  vector  of  the  blades  to  tilt 
forward  and  produce  a  propulsive  force  in  the  direction  of  rotation.  While  flapping 
the  blades  at  zero  collective  can  obviously  supply  a  rotational  force,  the  group  at 
Delft  has  shown  that  blade  flapping  eliminates  the  need  for  torque  into  a  helicopter 
rotor  at  positive  collectives  as  well.  At  a  blade  element  level,  the  effect  of  zero 
collective  and  positive  collective  is  shown  in  Figure  18  (a)  and  (b),  respectively.  The 
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first  case  (a)  produces  zero  net  lift  but  a  large  forward  thrust.  This  is  simply  the 
Knoller-Betz  effect.  The  second  (b)  demonstrates  positive  net  lift  and  positive  thrust, 
where  average  thrust  equals  the  average  drag.  It  is  in  this  regime  where  there  is  no 
torque  required  into  the  rotor  shaft.  The  “Omicopter”  may  be  the  only  flapping  rotor 
previously  tested  in  a  scientific  environment.  Because  of  the  large  structural  forces  of 
helicopter  blades  flapping  at  high  amplitudes,  this  idea  may  be  better  applied  to  an 
MAY. 


Figure  18:  Forces  on  flapping  blade 
at:  a)  zero  collective  and  b)  positive  collective 


A  micro  rotor  with  blades  that  can  independently  flap,  pitch,  and  rotate  was 
developed  and  patented  by  Bohorquez  and  Pines  [43],  but  it  was  not  tested  with  high 
frequency  flapping.  It  was  tested  in  pitching  to  a  frequency  of  two  flaps  per  rotor 
revolution.  It  was  found  that  higher  maximum  thrust  coefficients  and  slightly 
reduced  power  coefficients  were  possible  with  active  high  frequency  blade  pitching. 
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1.5  Current  Research  Objectives 


Methods  developed  in  nature  for  high  levels  of  lift  and  thrust  at  low  Reynolds 
numbers  can  also  be  applied  to  MAVs.  Biologically  inspired  methods,  such  as 
flapping  wings,  can  increase  the  thrust  of  MAVs  while  potentially  reducing  their 
power  required  to  fly.  An  increase  in  thrust  and  a  simultaneous  reduction  in  power 
would  lead  to  an  increase  in  payload,  a  more  compact  design,  and  an  increase  in 
endurance. 

The  objective  of  the  current  research  project  is  to  utilize  a  combination  of 
conventional  rotation  plus  avian  based  flapping  motion  to  increase  the  performance  of 
an  MAV  rotor.  Through  avian  based  blade  flapping  in  conjunction  with  passively 
flexible  blades,  this  flapping  rotor  should  be  able  to  reduce  or  even  eliminate  the 
torque  necessary  to  rotate  its  blades  and  produce  thrust.  This  MAV  flapping  rotor 
should  harness  delayed  stall  and  rotational  lift.  These  methods  should  be  introduced 
one  at  a  time  into  a  well  documented  system  in  order  to  understand  and  measure  their 
individual  effects.  This  dissertation  will  describe  in  detail  the  design  and  testing  of  a 
flapping  rotor  system  for  an  MAV. 

1.6  Outline 

To  fully  understand  the  problems  and  methods,  this  thesis  is  organized  as 
follows.  In  the  first  chapter,  performance  goals  for  an  MAV  are  defined.  An 
overview  of  state  of  the  art  technology  in  this  area  is  presented,  along  with 
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advantages  and  disadvantages  of  each  configuration.  Prior  related  research  and 
relevant  results  have  been  used  to  describe  aerodynamic  challenges.  Some  potential 
methodologies  by  other  researchers  for  overcoming  these  challenges  are  presented. 

Chapter  two  describes  a  proposed  solution,  which  is  a  small  flapping  rotor 
system  dubbed  the  “Flotor”.  The  advantages  and  disadvantages  of  such  an  idea  are 
discussed.  The  design  process  is  described  step  by  step  and  compared  to  empirical 
trends  from  nature.  The  initial  sizing  and  motion  requirements  are  presented, 
followed  by  the  design  and  fabrication  of  the  actual  research  test  bed. 

In  chapter  three,  the  design  is  analyzed.  The  geometry  and  kinematics  of  the 
mechanism  are  described  graphically  and  mathematically.  The  aerodynamic  analysis 
is  developed  using  quasi-steady  approximations  and  theoretical  calculations  are 
discussed. 

Chapter  four  presents  the  experimental  techniques  and  data  gathered  with  the 
test  bed.  The  experimental  parameters,  constraints,  and  assumptions  are  described 
and  the  results  are  presented. 

In  chapter  five,  conclusions  are  drawn  from  the  results  and  apparent  anomalies 
are  explained.  This  final  chapter  concludes  the  development  of  the  Flotor  flapping 
rotor  concept  for  micro  air  vehicles.  The  future  of  this  research  project  is  discussed  in 
terms  of  improvements  and  the  next  steps  for  demonstrating  the  concept  viability  for 
a  flying  vehicle. 
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2  Concept  and  Design 


Figure  19:  Flotor 


The  Flotor,  seen  in  Figure  19,  is  a  micro  rotor  test  stand  that  is  used  to 
demonstrate  the  effects  of  high  frequency  blade  flapping  on  average  rotor  forces.  It  is 
capable  of  independent  flapping  and  rotating  motion  and  it  can  rotate  freely  due  to  the 
aerodynamic  forces  caused  by  blade  flapping. 

2.1  Advantages 

The  Flotor  concept  combines  the  advantages  of  an  avian  based  flapper  and  a 
conventional  rotor.  Like  an  insect,  the  maximum  airspeed  of  the  blades  is  at  the  blade 
tips  and  the  blades  flap  about  a  root  with  near  zero  airspeed.  The  Flotor  could 
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produce  higher  thrust  coefficients  than  a  conventional  rotor  due  to  a  combination  of 
unsteady  and  conventional  aerodynamic  mechanisms,  delayed  stall  and  rotational  lift. 
Like  an  avian  based  flapper,  the  blades  of  the  Flotor  could  generate  their  own 
propulsive  force  in  the  chordwise  direction,  as  in  Figure  18.  The  propulsive  force 
from  flapping  will  cause  rotation  about  a  central  hub.  The  Flotor  can  utilize  the 
freestream  velocity  from  rotation  to  generate  lift,  like  a  bird  in  forward  flight.  To 
maintain  a  vertical  force,  the  blade  angles  on  the  Flotor  do  not  need  to  be  as  carefully 
trimmed  as  on  an  insect  or  avian  flapper.  The  benefit  of  constant  blade  rotation  about 
a  central  hub  is  that  the  aerodynamic  forces  are  easily  averaged  into  a  vertical  thrust 
vector,  resulting  in  a  hover  capable  device. 

2.2  Disadvantages 

Some  disadvantages  of  the  Flotor  are  also  similar  to  the  disadvantages  of 
flapping  wing  vehicles.  The  mechanism  for  blade  flapping  is  more  complex  than  a 
conventional  rotor  and  the  flapping  efficiency  is  not  well  established.  There  are 
periodic  inertial  forces  that  are  directly  proportional  to  the  mass  of  the  blades.  Since 
the  blades  in  this  case  flap  vertically,  the  inertial  forces  oscillate  in  a  vertical 
direction.  Aerodynamically  and  structurally,  the  blades  must  be  more  carefully 
optimized  than  blades  on  a  conventional  rotor. 

2.3  Initial  Rotor  Design  and  Sizing 

The  basic  blade  length  was  constrained  to  be  on  the  order  of  three  inches  in 
order  to  fit  within  the  accepted  size  range  for  an  MAY.  The  rotor  diameter  should  be 
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on  the  order  of  six  inches.  Starting  with  the  size  constraint,  the  remainder  of  the 
initial  design  was  compared  to  similar  size  flying  creatures.  Birds,  insects,  and  bats 
all  fit  tight  trends  of  sizing  and  performance  parameters  such  as  mass,  wing  span, 
aspect  ratio,  wing  loading,  flight  speed,  and  flapping  frequency.  The  parameter 
relations  can  usually  be  fitted  to  an  equation  of  the  form  y  =  amP  [44,45],  where  y  is 
the  parameter  in  question,  m  is  the  mass  of  the  flying  creature,  and  a  and  /3  are 
constants  fitted  to  the  data  with  a  least  squares  fit. 

2.3.1  Thrust 

In  nature,  creatures  with  a  wingspan  of  six  inches  are  generally  birds  or  bats, 
so  the  thrust  of  the  Flotor  can  be  compared  to  the  thrust  of  birds  and  bats  of  similar 
size.  Figure  20  [44]  shows  the  correlation  between  mass  and  wingspan  for  flying 
creatures  in  nature  as  well  as  the  DARPA  goal  for  an  MAV.  While  the  MAV  goal  is 
15  cm  span  and  100  grams  of  vehicle  mass,  flying  creatures  with  a  span  of  15  cm 
have  a  mass  of  10  grams  at  the  most.  The  relationship  between  the  mass  ( m )  and 
wingspan  ( b )  of  small  bats  fits  the  equation  b  =  1.3 m034,  which  corresponds  to  a  mass 
of  about  4  grams  for  a  span  of  6  inches  [44],  According  to  a  compilation  of  sizing 
parameters  for  all  birds  (excluding  hummingbirds)  [46,47,48],  the  relationship 
emerges  of  b  =  1.17m0  39.  This  corresponds  to  a  bird  mass  of  8  grams  for  a  span  of  6 
inches.  The  sizing  parameters  for  hummingbirds  show  the  relationship  2b  = 
2.24m0  53,  also  corresponding  to  a  mass  of  about  8  grams  [49]. 
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Figure  20:  Wingspan  vs.  body  mass  for  various  birds  and  bats  [44] 
and  sizing  goal  for  MAVs. 


Some  flying  creatures  can  produce  thrust  considerably  greater  than  their 
weight.  Some  bats,  for  example,  can  carry  their  young  in  flight  and  experiments  have 
shown  that  they  will  willingly  take  off  and  fly  even  when  loaded  with  over  70%  of 
their  bodyweight  on  their  backs  [50],  In  constant  forward  flight,  they  would  be  able 
to  carry  considerably  more.  For  such  reasons,  the  Flotor  should  be  able  to  produce 
more  thrust  than  just  the  weight  of  similarly  sized  flying  creatures.  The  high  thrust 
goals  compared  to  biological  flyers  requires  a  combination  of  higher  wing  loadings, 
higher  lift  coefficients,  and  higher  flight  speeds,  therefore  requiring  more  power  and 
possibly  reduced  endurance  [22], 
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2.3.2  Blades 


The  first  blades  manufactured  for  the  Flotor  were  sized  to  imitate  the  wings  of 
flying  creatures  of  similar  size.  Their  shape  and  scale  are  similar  to  the  wings  of 
small  birds  or  bats  and  their  structure  is  similar  to  that  of  a  bat.  Based  on  mass,  the 
aspect  ratio  relationship  for  bats  is  b2/S  =  11m011  [44]  and  for  birds  is  b2/S  = 
8.56 m0  06  [46],  According  to  this  relationship,  bats  of  this  size  have  wings  with  an 
aspect  ratio  of  about  6.5,  and  birds  have  wings  with  an  aspect  ratio  of  about  6.6.  In 
nature,  as  the  mass  decreases,  so  does  the  aspect  ratio. 

Since  the  total  blade  span,  the  aspect  ratio,  and  the  thrust  are  relatively 
constrained,  the  blade  area  and  blade  loading  can  be  directly  calculated  and  also 
compared  to  trends  in  nature.  With  a  blade  span  of  16  cm  and  an  aspect  ratio  of  6.5, 
the  blade  area  will  be  40  cm2.  With  a  thrust  of  about  20  grams  force,  the  wing 
loading  is  0.5g/cm2.  This  level  of  wing  loading  is  about  two  times  higher  than  that  of 
flying  creatures  of  this  size  [48],  In  general,  larger  flying  creatures  attain  higher  wing 
loadings,  but  they  also  fly  at  higher  airspeeds. 

2.3.3  Rotational  Speed 

In  general,  avian  flight  requires  a  forward  airspeed  for  the  wings  to  generate 
lift.  Except  for  hummingbirds,  birds  and  bats  utilize  forward  flight  for  their  most 
efficient  flight.  The  empirical  models  for  the  small  bird  group  “Passeriforms”  show 
that  the  speed  for  maximum  range  follows  the  trend  8.63 m0M.  Bats  flying  long 
distances  tend  to  fly  at  a  velocity  based  on  their  mass  of  1 1.8m  '  [45].  For  a  flying 

creature  in  the  size  regime  with  the  Flotor,  these  empirical  relations  point  to  a  forward 
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flight  speed  (V)  of  5  or  5.2  m/s,  respectively.  For  comparison,  the  reference  airspeed 
of  the  Flotor  blades  is  calculated  at  3A  of  the  rotor  radius.  Since  the  Flotor  radius  (R) 
is  0.08m,  the  Flotor  should  rotate  its  blades  at  approximately 
v  5  lm/s 

Q  =  —  =  — - =  85rad/s ,  or  81 1  RPM. 

\ R  0.06m 

2.4  Flapping 

The  flapping  frequency  and  flapping  amplitude  of  biological  flyers  generally 
decreases  with  increased  size.  Relationships  comparing  flapping  frequency  with 
body  mass  are  seen  in  Figure  21  and  show  a  great  deal  of  scatter  [51,52],  The 
predicted  upper  and  lower  flapping  frequency  limits  for  a  creature  with  a  mass  on  the 
order  of  10  to  20  grams  is  in  the  range  of  10  to  50  Hz.  Since  this  parameter  has  such 
a  wide  range  of  variation  in  nature,  it  will  be  one  of  the  main  variables  in  the  analysis 
and  experiment. 


Body  mass  (kg) 

Figure  21:  Wingbeat  frequency  vs.  Body  mass  for  birds  [51,52] 


Flapping  provides  propulsion  in  forward  flight,  but,  at  higher  frequencies,  it 
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also  introduces  unsteadiness  into  the  airflow.  The  measure  of  this  unsteadiness  is 


commonly  referred  to  as  “reduced  frequency”,  k,  where  k  =  .  The  flapping 

frequency  is  co,  the  wing  chord  is  c  and  the  airspeed  is  U.  The  flow  is  assumed  steady 
when  /c<0.05.  When  k> 0.2,  the  flow  is  so  unsteady  that  the  time  varying  airloads  can 
no  longer  be  accurately  predicted  by  conventional  steady  or  quasi-steady 
aerodynamics  [53].  Biological  flyers  in  the  size  range  of  the  flotor  operate 
completely  in  the  unsteady  regime.  The  general  reduced  frequency  trend  in  nature 
can  be  estimated  with  the  parameters  that  have  already  been  provided.  For  a  creature 
of  the  size  of  the  Flotor  with  a  wingspan  of  0.16  m  and  a  total  area  of  0.004  m2,  the 
average  wing  chord  is  about  0.025m.  The  forward  flight  speed  is  5.1m/s.  The  range 
of  flapping  frequencies  is  between  10  and  50  Hz,  or  63  and  314  rad/s.  These 
parameters  mean  the  reduced  frequency  for  flying  creatures  in  this  regime  generally 
lies  between  0.15  and  0.77,  excluding  hummingbirds.  Some  hummingbirds  flap  their 
wings  even  faster.  Thrush  Nightingales,  for  example,  are  a  bird  with  a  mass  of  about 
30  grams.  They  fly  with  a  reduced  frequency  of  0.19  when  in  fast  forward  flight  and 
0.53  when  in  slow  flight  [54],  The  Flotor  is  designed  to  operate  in  this  range. 

2.5  Design  and  Fabrication 

The  Flotor  was  designed  as  an  experimental  test  bed  capable  of  independently 
flapping  and  rotating  its  blades.  The  mechanism  was  designed  to  be  a  robust  platform 
with  predictable  movement  and  easily  adjustable  motions.  It  was  configured  in  order 
to  accurately  measure  the  aerodynamic  forces  created  during  experiments.  It  was  not 
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designed  for  efficiency  of  flapping  or  minimum  weight.  Because  low  weight  was  not 
a  priority,  the  safety  factor  of  most  parts  is  an  order  of  magnitude  higher  than  it  would 
be  for  a  flying  vehicle.  Most  of  the  Flotor  was  constructed  with  aluminum  due  to  its 
easy  machinability  and  relatively  high  strength  and  rigidity.  Off  the  shelf 
components  were  utilized  when  possible  to  reduce  the  time  for  manufacture. 

2.5.1  Basic  Components 

The  blade  flapping  actuation  of  the  Flotor  is  based  on  blade  pitch  actuation 
in  full  scale  helicopters.  A  ball  bearing  swashplate  is  utilized  to  transfer  loads  from 
the  fixed  frame  to  the  rotating  frame.  Instead  of  actuating  the  blade  pitch,  the 
swashplate  directly  actuates  the  blade  flapping  collectively.  Moving  the  swashplate 
vertically  causes  the  blades  to  flap  in  phase  in  the  same  direction.  The  rotor  and 
swashplate  assembly  is  shown  in  Figure  22. 

The  links  from  the  swashplate  to  the  blade  grips  are  attached  along  the  blade 
chordwise  center  of  mass  so  that  the  periodic  forces  from  the  links  do  not  cause  any 
flap-feather  coupling.  The  swashplate  link  ends,  which  are  hidden  inside  the  blade 
grips,  contain  ball  bearings  to  sustain  the  flapping  load  with  minimum  friction.  The 
lower  end  of  the  links  are  attached  to  the  swashplate  with  a  pin,  since  that  connection 
rotates  only  a  small  amount.  The  swashplate  and  rotor  shaft  operate  concentrically 
and  without  contact.  The  swashplate  is  constrained  to  linear  vertical  motion  with  a 
low  friction  Tusk  BX2-1  linear  slide. 
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Figure  22:  Rotor  and  swashplate  components 


A  two  view  diagram  of  the  linear  actuation  design  is  shown  in  Figure  23.  The 
non-rotating  swashplate  is  fastened  onto  the  linear  slide  with  screws,  making  a  firm 
cantilever  connection.  Also  attached  to  the  linear  slide  is  a  bracket  that  secures  a 
connecting  rod  with  ball  bearings,  making  a  pin  connection.  The  other  end  of  the 
connecting  rod  is  also  a  ball  bearing  pin  connection.  It  rotates  about  an  offset  pin  on 
the  flywheel  of  a  Hacker  brushless  motor,  which  is  called  the  flapping  motor.  The 
mechanism  allows  a  rotary  motor  to  produce  linear  actuation  with  very  little  friction. 

When  the  swashplate  is  actuated,  the  flapping  motion  of  the  blades  is  allowed 
by  a  spring  steel  flexbeam.  The  flexbeam  is  sandwiched  by  the  two-piece  hub  in  the 
center  and  by  the  blade  grips  on  both  ends.  The  clamping  force  for  the  hub  and  the 
blade  grips  comes  from  two  screws  that  pass  through  the  flexbeam  at  each  attachment 
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point.  With  a  blade  length  of  3  inches,  this  configuration  causes  a  hinge  offset  of  5% 
and  a  root  cutout  of  25%,  as  seen  in  the  drawing  in  Figure  24. 


The  shaft  of  the  rotor  has  two  ball  bearings  at  its  base,  which  constrain  its  axis 
of  rotation.  The  rotor  can  be  operated  with  or  without  a  motor  to  power  the  shaft 
rotation.  Without  a  motor  in  contact  with  the  shaft,  the  rotor  can  freely  rotate, 
experiencing  aerodynamic  drag  and  minimal  ball  bearing  friction.  In  the  freely 
rotating  configuration,  aerodynamic  forces  from  the  flapping  blades  are  the  only 
forces  that  can  cause  the  rotor  to  rotate.  To  better  control  the  rotation  speed  of  the 
shaft,  a  Firefly  coreless  DC  motor  with  4: 1  planetary  gearhead  is  used  to  power  the 
shaft  rotation  in  most  cases.  This  is  called  the  rotation  motor.  The  rotation  motor 
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gearhead  output  shaft  drives  the  rotor  shaft  with  a  1 : 1  ratio  gear.  The  1 : 1  gear  ratio 
was  limited  by  size  constraints  of  the  design. 


root  cutout 


Figure  24:  Front  and  Top  view  of  swashplate  and  rotor 


2.5.2  3-D  Solid  Model  and  Construction 

Before  a  single  piece  was  constructed,  the  entire  flotor  went  through  multiple 
iterations  of  3-D  solid  modeling.  CATIA  was  chosen  for  its  ease  of  operation  and 
accessibility.  Each  part  was  designed  along  with  its  neighboring  parts  and  assembled 
in  3-D  computer  workspace  to  ensure  an  exact  fit.  The  parts  were  then  machined  in- 
house  using  a  miniature  lathe  and  a  three  axis  milling  machine,  both  with  CNC 
capability.  The  CNC  was  primarily  utilized  on  the  mill  for  accurately  cutting  2-D 
profiles  to  different  depths.  Figure  22  shows  an  iteration  of  part  of  the  3-D  model. 
The  completed  3-D  model  is  shown  in  Figure  25.  The  cylinder  on  the  rear  of  the 
flapping  motor  is  a  device  to  measure  torque.  A  catalog  of  parts  is  kept  for  posterity 
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in  the  Alfred  Gessow  Rotorcraft  Center  graduate  student  office. 


Figure  25:  3-D  model  of  Flotor 


3  Analytical  Models 

Prior  to  its  operation,  the  kinematics  of  the  Flotor  were  planned  and  modeled. 
The  kinematic  model  was  used  to  predict  the  blade  motion,  which  was  used  to  predict 
aerodynamic  forces.  From  the  forces,  the  overall  thrust,  power,  and  figure  of  merit 
can  be  predicted. 

3.1  Kinematics  and  Mechanics 

The  kinematics  of  the  Flotor  can  be  analyzed  by  breaking  it  down  into  more 
abridged  components.  First,  the  flapping  mechanism  will  be  analyzed,  starting  with 
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the  flapping  motor  and  flywheel  with  offset  pin.  The  mechanism  has  been  simplified 
in  Figure  26.  The  swashplate  has  not  been  modeled  as  a  separate  part,  because,  for 
the  purpose  of  this  simplified  model,  the  swashplate  is  rigidly  attached  to  the  linear 
slide.  This  simplification  is  for  predicting  the  flapping  motion,  so  the  swashplate 
rotation  is  not  included  in  this  diagram. 
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The  flapping  motor  flywheel  rotates  at  an  angular  velocity  a>.  The  connecting 
rod  of  length  L  is  pinned  to  the  flywheel  at  a  distance  C  from  the  flywheel  center. 
The  connecting  rod  has  a  pin-pin  connection  to  the  flywheel  pin  and  to  the  linear 
slide,  so  the  force  through  the  connecting  rod  is  only  in  tension  or  compression.  The 
force  through  the  connecting  rod  drives  the  linear  slide,  which  rolls  linearly  on  ball 
bearings.  The  position,  X,  of  the  linear  slide  above  the  flywheel  center  depends  on 
the  angular  position  of  the  flywheel,  8. 

The  swashplate  link  is  also  a  pin-pin  connection,  with  one  end  pinned  to  the 
swashplate  and  the  other  end  pined  to  the  blade  grip.  Because  of  very  small  angular 
deflections  of  the  swashplate  link,  the  vertical  distance  from  the  linear  slide  to  the  pin 
connection  at  the  blade  grip  can  be  assumed  constant.  When  the  linear  slide  moves 
vertically,  the  swashplate  link  is  assumed  to  move  vertically  by  the  same  amount  and 
carry  the  forces.  The  axial  force  through  the  swashplate  link  causes  a  moment  in  the 
flexbeam.  This  moment  bends  the  flexbeam  and  causes  the  rigidly  modeled  blade 
and  blade  grip  to  deflect  at  an  angle  p  with  respect  to  horizontal.  The  motion  of  these 
parts  will  now  be  described  and  simplified. 

First,  the  vertical  swashplate  motion  will  be  described.  The  linear  mechanism 
that  drives  the  swashplate  has  motion  much  like  a  piston  in  an  internal  combustion 
engine  and  actually  has  the  same  equations  of  motion  [55].  The  equation  for 
swashplate  position  is  solved  through  simple  trigonometry.  The  vertical  velocity  and 
acceleration  of  the  swashplate  is  then  found  by  differentiating  the  position  equation 
with  respect  to  time. 
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Displacement  X  =  Ceos 8  +  VT2  -  C2  sin2  8 


Acceleration 


X  = 


Ceos  8 


From  the  equations  alone,  it  can  be  seen  that  the  vertical  motion  of  the 
swashplate  is  not  sinusoidal.  In  the  case  of  the  Flotor,  the  connecting  rod  length  L  is 
2  inches  and  the  crank  radius  C  is  variable.  The  ratio  of  C/L  is  referred  to  as  the 


“crank  ratio”.  One  cycle  of  swashplate  motion  is  plotted  in  the  following  figures  for 
two  different  values  of  crank  ratio.  The  swashplate  displacement  is  shown  in  Figure 
27.  In  this  case,  the  motion  of  the  swashplate  for  a  crank  ratio  of  0.1  is  nearly 
indistinguishable  from  sinusoidal  motion.  The  larger  crank  ratio  of  0.5,  however, 
causes  an  obvious  variation  from  sinusoidal  motion  that  appears  to  be  a  shift  in  phase. 

The  variation  from  sinusoidal  motion  becomes  even  more  obvious  when 


examining  the  vertical  velocity  and  acceleration  of  the  swashplate.  The  velocity  is 
shown  in  Figure  28  and  the  acceleration  is  shown  in  Figure  29  for  two  different  crank 
ratios.  This  type  of  movement  will  be  apparent  in  any  mechanism  that  relies  on  a 
crank  and  connecting  rod  to  actuate  a  linear  motion.  It  is  important  to  take  into 
account  this  motion  when  analyzing  MAV  flapping  mechanisms,  which  sometimes 
utilize  large  crank  ratios  in  the  interest  of  compactness. 
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Figure  27:  Swashplate  displacement 
for  crank  ratios  of  a)  50%  and  b)  10% 


Figure  28:  Swashplate  vertical  velocity 
for  crank  ratios  of  a)  50%  and  b)  10% 
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Since  the  Flotor  mechanism  was  not  designed  with  a  strict  goal  of 
miniaturization,  it  can  utilize  a  small  crank  ratio  with  a  relatively  long  connecting  rod. 
The  crank  distance  C  on  the  Flotor  is  never  more  than  10%  of  the  connecting  rod 
length  L,  so  C/L=0.1  demonstrates  a  worst  case  scenario  for  the  Flotor’s  variation 
from  sinusoidal  motion.  Except  for  very  large  amplitudes,  the  flapping  motion  of  the 
Flotor  can  be  modeled  as  sinusoidal. 


Figure  29:  Swashplate  vertical  acceleration 
for  crank  ratios  of  a)  50%  and  b)  10% 


3.2  Aerodynamics 

The  aerodynamic  analysis  of  the  Flotor  uses  conventional  helicopter  theories 
with  some  modifications  for  low  Reynolds  number  effects  and  high  amplitude 
flapping.  The  first  analytical  model  was  for  a  conventional  MAV  rotor.  The  second 
model  introduces  flapping  into  the  quasi-steady  model. 
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3.2.1  Rotation 


In  a  blade  element  analysis,  the  blade  is  broken  into  sections  and  each  section 
is  analyzed  individually.  The  aerodynamic  forces  from  each  element  are  then 
summed  to  obtain  the  overall  forces  on  the  rotor.  This  analysis  is  performed  at  each 
collective  angle  from  0°  to  36°  in  steps  of  2°.  The  basic  radial  distances  and  in  plane 
velocities  on  a  helicopter  blade  element  are  shown  in  Figure  30. 


UR 


Figure  30:  Blade  and  blade  element  distances  and  velocities  [53] 


The  overall  rotor  radius  is  R  with  a  root  cutout  of  r#.  The  blade  is  rotating 
with  an  angular  velocity  of  Q  radians  per  second.  The  blade  element  width  is  dy  and 
is  located  a  distance  y  from  the  rotor  center.  The  in-plane  wind  velocity  due  to 
rotation  is  Q y.  In  hover,  it  is  assumed  that  the  radial  component  of  the  velocity,  Ur, 
is  zero.  A  cross  section  of  the  blade  element  and  the  aerodynamic  vectors  and  forces 
acting  upon  it  are  shown  in  Figure  3 1 . 
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Figure  31:  Blade  element  angles,  vectors,  and  forces 


Uj  is  the  airspeed  of  the  blade  section  in  the  plane  of  rotation.  Up  is  the 
average  inflow  velocity  normal  to  the  plane  of  rotation.  The  total  resultant  velocity 

experienced  by  the  blade  section  is  U  =  ^U p2  +  U T2  .  0\s  the  collective  angle  of  the 

blade  section.  ^  =  tan  is  the  inflow  angle  and  a=  6  -  (j)  is  the  blade  section 

angle  of  attack.  From  these  vectors  of  wind  velocities,  the  lift,  dL,  and  drag,  dD, 
vectors  on  each  blade  section  can  be  estimated  and  added  to  get  the  resultant  force 
dFR.  The  horizontal  and  vertical  forces,  dFx  and  dFz,  are  derived  from  breaking  the 
resultant  force  vector  into  two  orthogonal  components. 

For  full  scale  helicopter  blade  element  analysis,  three  simplifying  assumptions 
are  usually  made.  First,  the  inflow  velocity  Up  is  assumed  much  smaller  than  UT,  so 
that  U  is  approximated  equal  to  UT.  Second,  the  induced  inflow  angle  <j)  is  assumed 
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small,  so  that  </>~  Up/Ut,  sin  <j>~  <f>,  and  cos  <f>~  1.  Third,  the  drag  is  at  least  one  order 
of  magnitude  smaller  than  the  lift,  so  the  vertical  contribution  due  to  drag,  dDsinfi,  is 
negligible  compared  to  the  vertical  contribution  due  to  lift  [53].  In  order  to  achieve 
an  accurate  model  of  rotor  forces  from  blade  element  analysis,  an  assumption  must  be 
made  for  the  inflow  distribution.  The  initial  estimate  of  inflow  will  utilize  the 
previous  three  assumptions. 

The  blade  element  momentum  theory  combines  the  blade  element  theory  and 
the  momentum  theory  to  predict  inflow  and  load  for  each  blade  element.  This  type  of 
analysis  can  be  utilized  to  predict  an  inflow  distribution  for  any  rotor  geometry. 
During  the  blade  element  momentum  model  used  for  the  current  analysis,  the  rotor 
disk  is  broken  up  into  20  concentric  annuli,  such  as  the  one  shown  in  Figure  32,  and 
the  blades  are  analyzed  as  20  corresponding  elements. 


Figure  32:  Annulus  used  for  momentum  theory  analysis  [53] 
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The  thrust  on  each  element  and  annulus  is  described  separately  by  blade 
element  theory  and  by  momentum  theory.  The  rotor  has  Nb  number  of  blades  with  a 
chord  of  c.  In  blade  element  theory,  the  thrust,  dT,  and  thrust  coefficient,  dCT,  are 
based  on  the  equation  for  blade  element  lift,  dL.  The  blade  section  is  operating  at  a 

lift  coefficient  of  C/.  Using  the  simplification  r  =  —  =  and  incorporating 

R  OR 

N  c 

solidity  cr  =  — — ,  the  basic  equation  for  the  thrust  coefficient  in  the  blade  element 
ttR 

theory  is  derived: 

1  2 

dL  =  —  pU  cC,dy 

dT  =  NbdL  =  Nb\  pU2cC,dy 

dc  NhjpU^cqdy 
T  p(ttR2)(QR)2 

dCT  =joC,r2dr 

The  momentum  theory  solution  for  thrust  is  the  change  in  momentum  flow 
through  each  rotor  annulus.  The  change  in  momentum  flow  across  each  rotor  annulus 
is  the  mass  flow,  dm  ,  multiplied  by  the  induced  velocity,  v„  multiplied  by  two.  The 
factor  of  two  is  due  to  the  rotor  wake  contraction,  which  causes  the  axial  flow  in  the 
free  wake  to  be  two  times  the  axial  inflow  at  the  rotor  disk  [53,  56],  The  inflow  ratio 
is  A=  Vi  /(Qy),  where  v,  =  Up  in  a  steady  state  with  no  blade  flapping.  The  derivation 
of  the  thrust  coefficient  in  hover  proceeds  as  follows: 
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dm  =  Inpvydy 
dT  =  Ivybh 


dT  =  4 npvi  ydy 


4npvi  ydy 
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The  inflow  for  each  rotor  element  is  estimated  by  equating  the  two 
interpretations  for  dCj  from  the  blade  element  model  and  the  blade  element 
momentum  model.  With  the  identity  A  =  (f>r  and  recalling  a  =  9-<j>,  and  C,  =C,  a , 

the  two  equations  can  be  combined  into  a  quadratic  equation  and  solved  as  follows: 


1  ,  oC , 

dCT  =  —  oC,r2dr=  °( Or 2  -  Ar)dr 

dCT  =  4  A2  rdr 


4/9  r dr  = 


4  A2  = 


oC, 

la 

~Y 

oC, 
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{Or2  -  Ar)dr 
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{Or)--±{A) 


.  uCi  0T1 

.  L  /  n  \  in 


A2+— HV-— ^{0r)  =  0 

o  o 


Mr)  = 
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/I  +  ——6r  - 1 


oC, 


The  only  unknown  in  this  solution  is  the  lift  curve  slope,  Cl  .  In  a  full  scale 
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helicopter  rotor  analysis,  the  lift  curve  slope  is  5.7/rad.  For  Reynolds  numbers  on  the 
order  of  104,  however,  the  lift  curve  slope  is  generally  about  10%  less,  verified  by 
CFD  data  and  experimental  data  of  airfoils  and  wings  [16,  17,  18,  20],  Utilizing  the 
lift  curve  slope  to  estimate  the  inflow  is  a  method  that  is  only  valid  below  stall,  where 
the  lift  curve  is  linear.  Especially  when  the  angle  of  attack  of  the  blade  element 
reaches  into  the  nonlinear  region  of  the  lift  curve,  the  inflow  model  depends  on  the 
actual  element  vertical  force. 

The  vertical  and  horizontal  forces  on  each  blade  element  are  calculated  by 
finding  the  angle  of  attack  and  corresponding  lift  and  drag.  The  lift  and  drag 
coefficients  for  a  specific  airfoil  are  most  accurately  found  directly  from  a  table  of  the 
airfoil  data.  For  this  particular  hover  analysis,  it  only  makes  sense  to  choose  an 
airfoil  that  performs  well  at  this  Reynolds  number,  such  as  a  thin  cambered  plate. 
Lift  and  drag  data  from  Sunada  et.  al.  [18]  and  Okamoto  and  Azuma  [37]  for  thin 
cambered  plates  at  low  Reynolds  numbers  are  compliled  in  Figure  33.  The  separate 
data  both  describe  a  5%  thick  plate  cambered  9%.  The  test  labeled  “Azuma”  was 
completed  at  a  Reynolds  number  of  7600  in  a  wind  tunnel.  The  test  labeled  “Sunada” 
was  done  in  a  water  tunnel  at  a  Reynolds  number  of  4000.  The  solid  lines  are  a 
piecewise  fit,  utilized  by  the  present  rotor  analysis  so  that  the  lift  and  drag  are  both 
calculated  as  functions  of  angle  of  attack. 
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C|  &  Cd  vs  angle  of  attack  for  cambered  plate 


a  (deg) 

Figure  33:  Cambered  plate  force  coefficients  at  low  Reynolds  number  [18,  37] 


From  this  point  forward,  when  calculating  the  thrust,  drag,  and  inflow  of  the 
rotor  elements,  the  forces  and  velocities  are  resolved  into  vertical  and  horizontal 
forces  as  shown  in  Figure  31.  The  vertical  component  of  the  drag  is  taken  into 
account  in  the  thrust  calculation,  even  though  it  is  usually  neglected  in  full  scale 
predictions.  The  forces  on  each  blade  element  are  multiplied  by  the  number  of 
blades,  to  get  the  total  vertical  force,  dFz  and  the  total  horizontal  force,  dFx,  for  each 
rotor  section: 

dFz  =  Nb  ( dL  cos  (/>  -  dD  sin  <j>) 
dFx  =  Nb  {dL  sin  (f)-\-dD  cos  (/>) 

This  is  where  the  current  prediction  differs  from  conventional  rotor  analysis. 
Using  a  part  of  momentum  theory  usually  applied  to  the  analysis  of  a  propeller,  there 
is  also  an  angular  velocity  component  of  the  inflow,  vsw,  which  is  a  direct  result  of  the 
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torque  of  the  rotor.  The  torque  required  for  each  rotor  section,  dQ,  is  equivalent  to 
the  change  in  angular  momentum  through  the  rotor  annulus  multiplied  by  the  annulus 
radius.  The  angular  momentum  of  the  inflow  at  the  rotor  disk  is  half  of  its  value  in 
the  far  wake  downstream  of  the  rotor,  just  as  the  vertical  inflow  at  the  rotor  disk  is 
half  the  far  wake  value  [57,  58],  This  leads  to  the  momentum  theory  results  for  the 
swirl  velocity  and  the  torque.  The  torque  can  also  be  described  in  terms  of  the 
horizontal  force,  dFx,  multiplied  by  the  element  radius,  y.  The  equation  relating  dFx 
and  swirl  velocity  then  simplifies  to  the  same  form  as  the  relationship  between  the 
thrust  and  the  vertical  inflow: 

dQ  =  2vswydm  =  47rpvtvswy2dy 

dFx  =  4npvivswydy 


dFx  dFz 

4m’,vswydy  47rpv2ydy 

The  momentum  theory  determines  that  the  square  of  the  vertical  inflow 
velocity,  vt2,  is  directly  proportional  to  the  vertical  force  on  each  rotor  annulus,  dFy. 
Likewise,  the  swirl  velocity  multiplied  by  the  inflow  velocity,  vswvi  is  directly 
proportional  to  the  horizontal  force,  dFx.  The  swirl  velocity  can  then  be  calculated  in 
terms  of  the  inflow  velocity: 


v 


dFy 

dFz 


When  analyzing  full  scale  rotors,  the  swirl  velocity  is  easily  neglected  because 
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it  is  at  least  an  order  of  magnitude  less  than  all  other  wind  velocities.  For  MAV 
rotors  with  lower  tip  speeds,  lower  lift  to  drag  ratios,  and  higher  thrust  coefficients, 
the  swirl  velocity  can  perceptively  reduce  the  blade  element  velocity  due  to  rotation, 
Ut.  For  the  case  of  the  current  MAV  rotor  prediction,  Ut  is  directly  reduced  by  the 
swirl  velocity,  so  that: 


UT 


Q y  -  v, 


X 

dFz 


The  vertical  force  for  each  rotor  element  is  used  to  recalculate  the  vertical 


inflow  and  the  horizontal  force  is  used  to  calculate  the  swirl  inflow  based  upon  the 
momentum  theory  result.  The  inflow  velocities  are  then  used  to  recalculate  the  blade 
element  forces.  This  iterative  process  is  repeated  until  the  inflow  estimate  changes 
less  than  2%  from  one  iteration  to  the  next.  The  values  for  vertical  inflow,  swirl 
velocity,  thrust  and  torque  usually  converge  within  a  few  iterations. 

After  the  inflow  is  finalized  for  each  rotor  annulus,  a  tip  loss  factor  is 
introduced.  The  tip  loss  factor  accounts  for  some  of  the  3-D  effects  near  the  rotor 
tips,  which  are  dominated  by  tip  vortices.  This  correction  factor  has  the  effect  of 
increasing  the  inflow  velocity  near  the  blade  tips  without  increasing  the  rotor  thrust. 
The  equations  utilized  for  tip  loss  calculations  were  described  by  Prandtl  and  Betz 
[59],  later  modeled  by  Goldstein  [60]  and  Lock  [61]  and  further  simplified  by 
Leishman  [53]  for  helicopter  applications.  The  inflow  “correction  factor”,  F, 
decreases  as  the  blade  element  radial  position  increases.  F  is  a  function  of  the  inflow 
and  the  radial  position  of  each  blade  element.  Since  the  inflow  depends  on  F  and  vice 
versa,  the  calculation  is  carried  out  iteratively  as  follows: 
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The  inflow  is  derived  from  equating  the  thrust  from  blade  element  theory  to 
the  thrust  from  momentum  theory.  The  inflow  estimate  and  tip  loss  factor  are 
functions  of  the  lift  curve  slope,  which  means  this  estimate  is  less  accurate  at  high 
angles  of  attack.  Only  for  the  purpose  of  the  tip  loss  estimation,  the  inflow  is 
artificially  limited  to  a  maximum  value.  As  the  collective  increases,  if  the  thrust  of 
the  blades  is  no  longer  increasing,  the  tip  loss  does  not  increase.  The  series  of 
equations  converges  within  a  few  iterations.  The  final  step  is  adding  the  forces  of  the 
blade  elements  to  calculate  the  total  rotor  thrust,  torque,  and  figure  of  merit.  These 
calculations  are  done  for  each  collective  angle. 
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3.2.2  Flapping  and  Rotation 

For  adding  flapping  motion,  a  time  component,  t,  is  introduced  into  the  model. 
The  total  time  period  of  analysis  is  equivalent  to  one  blade  flapping  cycle,  broken  into 
100  steps.  The  blade  flapping  is  assumed  sinusoidal  at  a  frequency  of  cof  radians  per 
second  and  an  amplitude  /3max,  so  that  the  flapping  angle  is  /?  =  /?max  sin  coFt .  Each 

blade  section  has  a  vertical  position  and  velocity  due  to  flapping: 

hF  =  y  sin  /?max  sm(cort) 

VF  =  coFy  sin  /?max  copco Ft) 

The  vertical  velocity  is  also  called  the  heaving  velocity.  A  downward  heaving 
velocity  causes  an  upward  wind  velocity,  which  increases  the  blade  section  angle  of 
attack.  Under  the  conditions  seen  in  Figure  34,  the  resultant  force  on  the  blade 
section  tilts  in  a  forward  direction  and  the  drag  force  contributes  to  the  vertical  force. 
At  this  instant  of  downward  heaving,  the  blade  section  is  producing  a  propulsive  force 
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and  is  producing  a  vertical  force  that  is  greater  than  the  sectional  lift.  With  high 
frequency  blade  flapping,  the  heaving  motion  is  an  order  of  magnitude  larger  than  the 
inflow  velocity,  so  the  inflow  velocity  is  not  shown  separately  in  this  figure.  The 
inflow  velocity  is  not  ignored  in  the  calculations,  however.  The  vertical  velocity 
component  becomes  Up  =  v;.  +  VF .  This  change  propagates  throughout  the  model, 

i  Up 

starting  with  </>  =  tan  — — . 


Making  a  couple  of  interim  assumptions,  it  can  be  shown  that  the  freestream 
angle  variation  and  therefore  the  angle  of  attack  variation  is  approximately  constant 
across  the  span  of  the  blade.  The  maximum  and  minimum  freestream  angle  occurs 
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when  the  heaving  velocity  is  also  at  its  maximum  and  minimum.  The  heaving 
velocity  is  approximately  VF  =  coFy  sin/?max  cos(®Ft),  so  the  maximum  heaving 

velocity  is  VF  =  coFy  sin  /?max .  Ignoring  the  swirl  velocity,  the  forward  velocity  is 


approximately  Ut=  Qy.  The  maximum  freestream  angle  becomes: 
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Therefore,  the  approximate  maximum  freestream  angle  depends  only  on  the  flapping 
frequency,  the  flapping  amplitude,  and  the  hub  rotation  speed.  It  should  remain 
relatively  constant  across  the  span  of  the  blades. 

There  is  an  option  in  the  model  so  that  the  blades  can  be  fundamentally 
modeled  as  flexible.  A  “twist  factor”  is  introduced  that  alters  the  blade  element  angle 
of  attack  depending  on  the  blade  element  load.  Since  the  blade  element  forces 
depend  directly  on  the  angle  of  attack,  the  angle  of  attack  simply  becomes  a  fraction 
of  its  previous  value,  depending  on  the  twist  factor. 

The  blade  element  analysis  of  the  forces,  velocities,  and  angles  are  calculated 
the  same  as  the  conventional  rotation  case  except  for  the  inflow.  The  inflow  velocity, 
Vi,  for  each  blade  element  and  collective  angle  is  borrowed  from  the  non-flapping 
case  and  is  assumed  constant  with  time.  Since  the  forces  from  this  flapping  case  are 
time  varying,  they  are  averaged  to  compare  them  with  the  results  from  the 
conventional  rotation  case. 
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4  Experiment,  Results,  and  Discussion 


The  experiments  with  the  Flotor  exploited  three  types  of  independently 
controlled  motion.  The  first  motion  tested  was  powered  flapping  motion  of  flexible 
blades  with  a  freely  rotating  main  rotor  shaft.  The  second  motion  was  conventional 
powered  rotation,  which  could  be  compared  with  previously  published  data  from 
similar  experiments.  The  third  and  final  group  of  tests  combined  simultaneous 
powered  rotation  and  powered  flapping  motions  of  stiff  blades.  Because  of  the  many 
variables  and  configurations  possible  with  the  test  apparatus,  the  variables  were 
narrowed  by  applying  rational  constraints  to  hold  constant  as  many  variables  as 
possible.  Each  set  of  tests  is  presented  with  its  own  variables  and  constraints. 

4. 1  Powered  Flapping  With  Free  Hub  Rotation 

In  the  tests  utilizing  flapping  with  free  shaft  rotation,  the  variables  included 
blade  design,  flapping  frequency,  and  flapping  amplitude.  The  hub  rotation  speed 
was  purely  a  function  of  these  three  variables  as  it  trimmed  itself  for  a  constant  speed 
with  zero  rotor  torque.  The  first  tests  with  the  Flotor  were  conducted  with  the 
apparatus  in  the  configuration  as  seen  in  Figure  35.  The  only  motor  in  this 
configuration  was  the  flapping  motor,  which  was  a  Hacker  B20  26L  brushless  motor 
with  4:1  planetary  gearhead.  The  brushless  design  required  the  use  of  an  electronic 
speed  controller.  To  control  the  speed  of  the  flapping  motor,  the  speed  controller  was 
input  a  pulse  signal  from  a  GWS  adjustable  pulse  generator  and  servo  tester.  The 
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flapping  motor  speed  was  varied  manually  by  modulating  the  pulse  width  output  from 
the  GWS  pulse  generator.  The  amplitude  of  flapping  was  fixed  by  the  distance  of  an 
offset  pin  in  the  flywheel  of  the  flapping  motor.  For  this  test,  the  amplitude  of 
flapping  was  ±  17°,  which  does  not  include  blade  flexibility. 


Figure  35:  Flapping  rotor  with  freely  rotating  hub 

For  the  free  rotation  of  the  rotor  to  be  self-starting,  it  was  necessary  for  the 
blades  to  pitch  into  the  relative  wind  when  they  started  flapping.  Therefore,  the 
blades  had  to  be  flexible  with  the  elastic  axis  located  in  front  of  the  aerodynamic 
center  of  pressure.  There  are  four  basic  cases  that  can  be  used  to  demonstrate  the 
necessity  for  this  type  of  flexibility:  rigid  blades  at  zero  collective,  rigid  blades  at 
positive  collective,  flexible  blades  at  zero  collective,  and  flexible  blades  at  positive 
collective.  All  of  the  cases  show  the  forces  on  a  flat  plate  blade  element  assuming 
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the  rotor  is  not  rotating  and  the  only  motion  is  blade  flapping. 

The  first  two  cases  of  rigid  blade  flapping  with  zero  and  positive  collective  are 
shown  in  Figure  36.  The  actual  velocity  of  this  blade  element  due  to  flapping  is  VF> 
which  means  the  only  wind  velocity  is  -VF  in  a  vertical  direction.  The  blade  element 
is  shown  at  the  point  of  maximum  velocity  during  the  downstroke  and  upstroke  of 
flapping,  assuming  that  the  rest  of  the  blade  motion  is  symmetric.  At  zero  collective, 
the  angle  of  attack  on  the  blade  is  90°  during  the  downstroke  and  the  upstroke.  A  flat 
plate  at  90°  angle  of  attack  produces  only  a  drag  force.  The  drag  force  during  the 
downstroke  and  the  upstroke  is  equal  and  opposite  and  therefore  averages  to  zero.  In 
an  inviscid  fluid,  there  is  no  force  in  the  plane  of  rotation  during  the  upstroke  or 
downstroke  at  zero  collective. 


Rigid  Blades 


Downstroke 


Upstroke 

-VP 


Zero 

Collective  VF 


Positive 

Collective 


Figure  36:  Flapping  rigid  flat  plate  blades 
at  a)  0°  collective  and  b)  Positive  collective 
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At  positive  collective,  the  rigid  blade  has  an  angle  of  attack  less  than  90° 
during  the  downstroke,  producing  a  lift  force  in  the  plane  of  rotation.  The  lift  force 
on  the  upstroke  and  the  downstroke  is  equal  and  opposite,  however,  and  averages  to 
zero.  The  drag  forces  are  also  equal  and  opposite.  Therefore,  flapping  rigid  flat 
blades  at  any  collective  angle  should  not  produce  a  net  force  in  any  direction  and 
should  not  cause  the  rotor  to  start  rotating. 

The  example  of  blades  flexible  about  their  leading  edge  is  shown  in  Figure  37. 
These  blades  are  assumed  to  have  their  elastic  axis,  marked  “ea”,  at  their  leading 
edge,  which  ensures  that  aerodynamic  forces  cause  the  blade  to  pitch  into  the  relative 
wind  during  flapping.  With  zero  degrees  average  collective,  the  aerodynamic  force 
pitches  the  blade  so  that  it  has  an  equal  angle  of  attack  during  the  upstroke  and  the 
downstroke.  The  drag  force  during  the  upstroke  and  downstroke  is  equal  and 
opposite  and  therefore  averages  to  zero  net  vertical  force.  The  lift  force,  however,  is 
in  the  same  direction  during  the  downstroke  and  upstroke,  producing  a  net  force  in 
the  plane  of  rotation.  This  net  force  accelerates  the  rotation  of  the  rotor.  At  a  positive 
collective  angle  with  flapping  flexible  blades,  the  angle  of  attack  on  the  downstroke  is 
always  greater  than  the  angle  of  attack  on  the  upstroke.  In  the  figure,  the  downstroke 
is  shown  with  nearly  90°  angle  of  attack,  producing  no  lift  in  the  plane  of  rotation,  but 
producing  a  high  vertical  drag  force.  During  the  upstroke,  the  blade  pitches  to  a 
much  more  acute  angle  of  attack,  producing  a  high  lift  in  the  plane  of  rotation  and 
relatively  little  drag  in  the  vertical  direction.  In  this  case,  the  lift  averages  to  a  net 
force  in  the  plane  of  rotation  while  the  drag  averages  to  a  net  vertical  thrust. 
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Figure  37:  Flapping  flat  plate  blades  flexible  about  their  leading  edge 
at  a)  0°  collective  and  b)  Positive  collective 


To  ensure  the  forward  location  of  the  elastic  axis,  the  blades  for  the 
experiment  were  manufactured  with  a  leading  edge  spar  that  is  an  order  of  magnitude 
stiffer  than  the  rest  of  the  blade.  The  blades,  shown  in  Figure  38,  were  manufactured 
out  of  three  beams  of  pre-cured  unidirectional  carbon/epoxy  composite  bonded  to  a 
0.25  mm  thick  sheet  of  mylar.  The  carbon  beams  have  fibers  running  in  their 
lengthwise  direction.  The  leading  edge  spar  is  the  stiffest  beam,  with  a  length  of  8.35 
cm,  thickness  of  0.65  mm  and  a  width  of  2.95  mm.  The  root  end  of  the  leading  edge 
spar  is  attached  to  the  hub  by  securing  0.45  cm  of  material  into  the  blade  grips.  The 
two  supporting  element  beams  are  more  flexible  due  to  their  reduced  dimensions  of 
0.25  mm  thick  and  2.9  mm  wide.  The  supporting  element  perpendicular  to  the 
leading  edge  spar  is  2.85  cm  long.  The  diagonal  supporting  element  is  positioned  at 
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an  angle  of  22°  from  the  leading  edge  and  extends  to  the  trailing  edge  at  60%  blade 


span. 


The  blades  used  for  this  test  have  a  lifting  surface  with  a  length  of  7.92  cm. 
Including  the  root  cutout,  which  is  30%  in  this  case,  the  rotor  radius  is  10.95  cm. 
There  are  two  linear  taper  ratios  incorporated  into  the  blades.  The  first  taper  ratio  is 
14%  and  extends  from  the  root  to  60%  span.  The  second  taper  ratio  of  55%  extends 
from  the  60%  span  to  the  tip,  which  gives  the  blades  a  shape  representative  of  a  bat  of 
similar  size  [44],  The  root  chord  of  a  blade  is  3.15  cm,  the  chord  at  60%  span  is  2.47 
cm,  and  the  chord  at  the  tip  is  0.63  cm.  The  total  area  of  both  blades  is  39.5  cm2,  or 
0.00395  m2,  which  does  not  include  the  root  cutout  of  the  rotor. 

The  purpose  of  the  first  test  was  only  to  demonstrate  the  concept  of  passive 
hub  rotation  with  powered  flapping.  The  variables  in  this  test  were  blade  collective 
angle  and  the  blade  flapping  frequency.  The  collective  angle  and  flapping  amplitude 
were  measured  with  a  bubble  type  inclinometer.  The  flapping  amplitude  was 
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measured  by  manually  positioning  the  flywheel  for  maximum  positive  and  negative 
deflection  and  taking  measurements  at  each  position.  The  flapping  amplitude  was 
fixed  at  ±17°,  or  34°  peak  to  peak.  In  the  first  test,  the  collective  was  varied  from  0° 
to  15°  to  30°. 

The  speed  of  the  flapping  motor  was  varied  between  26  and  43  Hz,  or  1600 
and  2600  RPM.  The  rotational  speeds  of  the  hub  and  flapping  motor  were  both 
measured  with  a  Monarch  Instruments  light  sensitive  tachometer.  Two  strips  of 
reflective  tape  were  placed  1 80°  apart  on  each  rotating  part.  The  tachometer  works 
by  emitting  a  beam  of  light  pointed  at  the  reflective  tape,  which  reflects  back  to  the 
tachometer.  The  tachometer  was  calibrated  to  count  two  pulses  of  light  for  each 
rotation  to  display  the  proper  rotational  speed. 

This  first  set  of  results  proved  that  the  hub  can  rotate  freely  due  to  the 
aerodynamic  forces  of  the  flapping  blades,  even  when  the  blades  are  at  high  collective 
angles.  The  results  gathered  during  this  proof  of  concept  test  were  preliminary  and 
are  shown  in  Figure  39.  At  0°  collective,  the  hub  rotation  varies  from  170  to  525 
RPM  as  the  flapping  varies  from  28  to  43  Hz.  At  15°  collective,  the  hub  rotation 
varies  from  190  to  500  RPM  as  the  flapping  speed  varies  from  27  to  43  Hz.  At  30° 
collective,  the  hub  rotation  varies  from  290  to  550  RPM  as  the  flapping  speed  varies 
from  27  to  43  Hz.  The  ratio  of  hub  rotation  speed  to  flapping  speed  is  not  constant. 
With  the  blades  at  15°  collective,  the  ratio  of  hub  rotation  to  flapping  speed  increases 
with  increased  flapping  speed  at  every  data  point.  At  0°  and  30°  collective,  the  ratio 
increases  at  all  but  the  highest  flapping  speed.  At  30°  collective,  the  hub  rotation 
speed  is  highest  for  every  flapping  frequency. 
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At  the  highest  rotational  speeds  of  550  RPM,  the  3A  chord  of  the  blades  is 
reaching  airspeeds  of  approximately  4.5  m/s.  Compared  to  similarly  sized  flying 
creatures,  the  blade  velocity  in  this  configuration  falls  short  of  its  target.  The  obvious 
trend  in  this  data  shows  that  increasing  flapping  frequency  is  related  to  increasing  the 
hub  rotation  speed.  The  flapping  frequency  in  this  case  was  limited  by  the  forces  of 
the  reciprocating  parts  on  the  measurement  equipment.  The  test  stand  was  mounted 
upon  a  thrust  cell  which  had  a  maximum  capacity  of  1  kg,  so  the  total  calculated 
inertial  forces  of  all  reciprocating  parts  was  kept  below  1  kg  force. 

The  first  average  data  from  the  thrust  cell  is  shown  in  Figure  40.  The  data 
was  gathered  to  ensure  all  components  of  the  Flotor  as  well  as  the  instrumentation  for 
data  acquisition  functioned  properly.  Although  preliminary,  these  results  showed 
increasing  thrust  with  increasing  flapping  frequency  at  every  data  point.  The  thrust 
increases  with  increasing  collective  up  to  30°,  but  beyond  30°  no  further  thrust 
increase  occurred. 
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In  order  to  measure  the  vertical  forces  from  the  rotor,  multiple  configurations 
to  gather  the  time  varying  thrust  data  were  attempted  but  none  were  successful.  Since 
the  inertial  loads  due  to  flapping  are  orders  of  magnitude  greater  than  the 
aerodynamic  forces,  it  is  a  very  difficult  task  to  differentiate  the  instantaneous 
aerodynamic  loads  from  inertial  forces.  However,  since  the  average  inertial  forces 
are  nearly  zero,  an  average  signal  will  contain  only  the  aerodynamic  force. 
Therefore,  only  the  average  forces  were  measured  during  the  remainder  of  the 
experiments.  To  measure  the  average  thrust,  the  test  stand  was  mounted  to  an  Ohaus 
Explorer  digital  balance.  The  balance  incorporates  a  feature  which  uses  a  low  pass 
filter  while  averaging  the  force  for  five  seconds  to  display  the  result.  It  eliminates  the 
periodic  flapping  forces.  This  feature,  which  is  normally  used  for  measuring  the 
weight  of  active  live  animals,  has  proven  very  useful  and  accurate  for  measuring  the 
average  aerodynamic  forces  of  the  Flotor. 
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Instead  of  a  manual  tachometer,  the  flapping  motor  and  rotor  hub  were 
instrumented  with  sensors  to  monitor  their  rotation.  The  RPM  data  was  gathered  with 
National  Instruments  hardware  and  software.  The  signals  from  the  sensors  were  fed 
through  SCC  modules  in  a  SC-2345  Carrier.  The  carrier  provided  power  to  the 
modules  and  also  interfaced  with  a  computer  through  a  6062E  data  acquisition  card. 
The  computer  ran  LabVIEW  software  for  data  acquisition.  This  setup  monitored 
inputs  and  provided  power  output.  The  power  output  in  this  case  was  also  used  to 
power  the  sensors  that  monitor  rotational  speeds. 

To  ensure  low  friction,  a  system  of  magnets  and  hall  sensors  was  devised  so 
that  no  contact  with  the  rotating  parts  is  necessary  to  measure  their  speed.  Because  of 
the  small  size  of  the  rotating  parts,  the  magnets  had  to  be  small  and  relatively 
powerful  to  activate  the  hall  sensors.  Neodymium  rare  earth  magnets  were  selected 
with  a  length  and  width  of  1/8”  and  a  depth  of  1/16”.  Two  of  the  magnets  were 
bonded  180°  apart  on  the  flapping  motor  flywheel  and  two  were  mounted  180°  apart 
on  the  rotor  hub.  Honeywell  SS443A  hall  effect  sensors  were  mounted  close  enough 
to  the  rotating  magnets  so  that  the  magnets  passed  within  a  few  millimeters  of  the 
face  of  the  sensors.  This  model  of  sensor  is  unidirectional,  so  it  was  important  that 
the  proper  pole  of  each  magnet  faced  the  sensor  as  it  passed.  The  sensor  was 
powered  by  a  constant  10  V  from  the  data  acquisition  module.  The  output  from  the 
sensor  is  a  “sinking”  output,  so  the  LabVIEW  software  was  written  to  count  one 
rotation  for  every  two  drops  in  voltage  from  the  sensor  output. 

The  next  test  was  the  first  that  simultaneously  incorporated  measurements  of 
hub  speed,  flapping  speed,  and  thrust.  The  rotor  was  operated  in  this  case  at  a  fixed 
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collective  of  15°,  which  previously  showed  relatively  high  thrust  levels  and  hub 
rotation  with  the  aforementioned  blades.  The  flapping  frequency  was  varied  between 
17  and  38  Hz. 

Similar  to  the  previous  test,  the  hub  RPM  again  increased  non-linearly  with 
blade  flapping  frequency.  Figure  41  shows  the  hub  rotating  speed  compared  to  the 
blade  flapping  frequency.  As  the  flapping  frequency  increased  from  16  to  38  Hz,  the 
rotor  hub  speed  nearly  quadrupled  from  150  RPM  to  550  RPM.  Figure  42  shows  the 
ratio  of  hub  to  flapping  speed  versus  the  flapping  frequency.  This  shows  that  as  the 
flapping  frequency  doubled,  the  ratio  of  hub  rotation  to  flapping  frequency  nearly 
doubled  as  well. 
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Figure  41:  Hub  RPM  vs.  flapping  frequency  for  15°  collective 

Figure  43  combines  the  thrust  and  hub  RPM  into  a  thrust  coefficient,  which  is 
plotted  versus  flapping  frequency.  The  thrust  coefficient  increases  by  40%  from  16  to 
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32  Hz  flapping  frequency.  At  higher  flapping  frequencies,  the  thrust  coefficient  Ct 
levels  off  and  even  decreases  slightly  at  the  highest  flapping  frequency.  These  thrust 
coefficients  are  on  the  order  of  two  times  as  much  as  the  thrust  coefficients  on  a  full 
scale  rotor. 
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Figure  42:  Hub  speed  ratio  vs.  flapping  frequency  for  15°  collective 


Increased  flapping  frequency  also  increases  the  rotational  speed  of  the  rotor 


hub.  Increasing  the  rotation  speed  of  a  conventional  rotor  at  a  certain  collective  has 


no  effect  on  the  thrust  coefficient.  In  this  case  at  15°  collective,  however,  the  thrust 


coefficient  increases  with  increasing  rotational  speed,  as  seen  in  Figure  44.  This 


increase  shows  that  the  blades  are  achieving  vertical  thrust  by  means  other  than 


steady  hub  rotation.  This  thrust  due  to  flapping  increases  the  total  thrust  as  flapping 


frequency  increases. 
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Figure  43:  Thrust  coefficient  vs.  flapping  frequency  at  15°  collective 
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Figure  44:  Thrust  coefficient  vs.  hub  RPM  at  15°  collective 


Because  of  the  low  tip  speeds  of  the  Flotor  and  most  MAVs,  the  average  blade 
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loading  of  the  Flotor  is  far  below  that  of  full  scale  helicopters.  The  highest  thrust 
during  this  test  was  about  4  g,  corresponding  to  a  blade  loading  of  about  10  N/m2, 
which  falls  exactly  in  the  range  of  bats  of  similar  size  [45], 

During  the  previous  test,  a  torque  cell  was  added  to  measure  the  instantaneous 
torque  of  the  flapping  motor.  The  rear  of  the  torque  cell  was  mounted  rigidly  and  the 
front  of  the  torque  cell  was  mounted  to  the  flapping  motor,  which  was  suspended  by 
ball  bearings.  Any  torque  produced  by  the  flapping  motor  was  transferred  through 
the  torque  cell,  then  recorded  by  the  data  acquisition  system.  The  torque  for  a  case 
with  32  Hz  flapping  is  shown  in  Figure  45.  For  each  flapping  cycle,  there  should  be 
two  torque  peaks,  one  while  accelerating  the  flapping  mechanism  upward  and  one 
while  accelerating  the  mechanism  downward.  The  torque  peak  during  the  upward 
acceleration  should  be  greatest,  since  it  is  opposite  in  direction  from  gravity.  This 
explains  the  largest  peaks  in  the  torque  signal.  The  aerodynamic  contribution  to  the 
torque  signal  at  this  frequency  and  amplitude  is  on  the  order  of  0.001  -  0.01  Nm.  This 
magnitude  is  barely  distinguishable  in  the  total  torque  reading.  In  any  case,  the 
aerodynamic  torque  should  peak  90°  after  the  mechanical  torque  because  the 
aerodynamic  forces  vary  with  flapping  velocity  while  the  mechanical  torque  varies 
with  flapping  acceleration. 
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Figure  45:  Torque  vs.  Time  for  32  Hz  Flapping 
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During  the  next  test,  the  hub  configuration  was  changed  slightly  to  reduce  friction. 
The  single  ball  bearing  that  was  previously  mounted  in  the  hub  was  replaced  with  two 
ball  bearings.  Even  so,  the  two  bearings  were  difficult  to  align  since  they  were 
pressed  into  holes  in  two  separate  parts.  After  this  test,  the  bearings  were  relocated  to 
the  base  of  the  shaft,  where  their  outer  rings  were  pressed  into  opposite  sides  of  a 
single  hole  and  the  shaft  pressed  through  the  center. 

The  flapping  frequency  was  varied  from  20-39  Hz  at  collective  angles  of  -20°, 
10°,  20°,  and  40°  and  the  shaft  was  free  to  rotate.  In  order  to  determine  the  thrusting 
of  the  blades  with  zero  rotation,  two  tests  were  done  with  the  shaft  rotation  locked. 
One  test  was  done  at  20°  and  another  at  40°  with  only  vertical  flapping  and  no  other 
motion.  In  these  tests,  the  blades  flapped  in  a  vertical  plane  with  the  same  range  of 
frequencies.  The  flapping  amplitude  was  34°  peak  to  peak. 

Each  data  point  shown  in  this  test  is  an  averaged  data  point  of  five  tests  of  the 
same  conditions.  Rotation  speed  data  from  the  test  with  blade  flapping  and  free  hub 
rotation  is  shown  in  Figure  46.  The  hub  rotates  the  fastest  at  10°  collective  at  all 
flapping  frequencies.  The  hub  rotation  with  the  blades  at  -20°  collective  is  very 
similar  to  the  hub  rotation  with  the  blades  set  at  20°  collective.  At  40°  collective,  the 
rotation  speed  is  the  slowest  at  every  flapping  frequency.  The  hub  rotation  again 
increases  at  a  greater  rate  than  the  flapping  speed  except  at  the  highest  flapping 
frequencies. 
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Figure  46:  Hub  rotation  vs.  flapping  frequency  at  different  collectives 


Figure  47  shows  the  thrust  coefficient  versus  flapping  frequency  for  different 
blade  collectives.  At  moderate  collective  angles,  the  thrust  coefficient  decreased  with 
flapping  frequency  because  of  its  inversely  proportional  relationship  to  the  square  of 
hub  rotation.  The  absolute  thrust  increased,  however,  due  to  the  increased  hub 
rotation  and  because  a  portion  of  the  thrust  was  produced  by  blade  flapping. 

— ♦ — 20  deg  — 1 0  deg  — 20  deg  x  40  deg 


Figure  47:  Thrust  coefficient  vs.  flapping  frequency  at  different  collectives 
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Even  though  the  thrust  coefficients  may  decrease  with  increased  flapping 
frequency,  the  absolute  thrust  values  increase  with  increased  flapping  frequency.  The 
absolute  values  of  thrust  are  shown  in  terms  of  blade  loading,  which  is  plotted  against 
flapping  frequency  in  Figure  48.  In  the  data  labeled  “no  rotation”,  the  thrust  is  due  to 
only  flapping  with  the  shaft  rotation  locked.  The  results  help  to  explain  why,  in  the 
previous  figure,  the  40°  collective  thrust  coefficient  increases  with  increased  flapping 
frequency.  The  thrust  from  flapping  without  rotation  at  40°  collective  is  up  to  40% 
of  the  thrust  of  flapping  with  free  rotation.  In  the  20°  collective  case,  however,  less 
than  20%  of  the  overall  thrust  can  be  attributed  to  the  blade  flapping.  The  majority  of 
the  thrust  in  the  moderate  collective  cases  comes  from  hub  rotation. 


Figure  48:  blade  loading  vs.  flapping  frequency 
with  rotation  and  with  no  rotation  at  20°  and  40°  collective 


The  results  in  all  of  the  previous  tests  maintained  clear  trends  at  all  but  the 
highest  flapping  frequencies.  When  the  flapping  frequency  reached  well  above  30 
Hz,  it  was  observed  that  the  blades  were  undergoing  extreme  deformations.  Apparent 
trends  that  appeared  at  lower  frequencies,  such  as  increased  hub  rotation  and 
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increased  thrust  with  increased  flapping  frequency,  changed  as  the  flapping  frequency 
reached  above  35  Hz.  It  was  thought  that  stiffer  blades  would  produce  a  more 
predictable  behavior  at  higher  flapping  frequencies. 

The  composition  of  the  next  blades  can  be  seen  in  Figure  49.  There  is  a 
leading  edge  spar  and  a  rectangular  blade  section  attached  to  the  spar  with  a  clear 
mylar  hinge.  The  material  chosen  for  the  leading  edge  spar  was  the  same 
unidirectional  carbon  epoxy  as  the  previous  blades.  The  leading  edge  hinge  behaves 
like  a  spring  during  flapping.  The  blade  material  itself  was  constructed  to  be 
relatively  stiff  compared  to  the  previous  blades  while  maintaining  light  weight.  The 
blade  surface  is  made  of  2  layers  of  carbon  epoxy  prepreg  weave  laid  at  ±45°  with 
respect  to  the  leading  edge.  This  gives  the  blade  fairly  stiff  properties  in  both  the 
spanwise  and  chordwise  direction  as  well  as  high  torsional  stiffness.  The  blade  has  a 
length  of  7.62  cm  and  a  chord  of  1.78  cm,  giving  the  two  blades  a  total  blade  area  of 
27.1  cm2,  or  0.0021  m2.  The  blades  give  the  rotor  a  total  radius  of  10.2  cm. 


Figure  49:  Carbon  blade  with  flexible  leading  edge  hinge 


These  blades  were  tested  at  collective  angles  ranging  from  -10°  to  30°  and  at 
flapping  frequencies  from  19  to  41  Hz.  The  collective  was  varied  in  steps  of  5° 
between  0°  and  30°  and  was  also  tested  at  -10°  to  verify  negative  thrust  readings.  The 
flapping  amplitude  was  the  same  as  the  earlier  tests  at  a  peak  to  peak  value  of  34°. 
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Like  the  previous  test,  each  data  point  shown  here  is  an  average  of  at  least  three  tests 
at  the  same  conditions.  Even  with  vastly  different  blade  constructions,  the  results 
were  similar  to  the  results  from  the  previous  blade  design. 

Figure  50  is  a  compilation  of  the  results  of  hub  rotation  speed  versus  flapping 
frequency  at  all  collective  angles.  Once  again,  the  blades  reached  a  maximum 
rotation  speed  at  15°  collective,  with  similar  results  at  20°  collective.  The  rotation 
speed  again  increased  at  a  faster  rate  than  the  flapping  frequency.  For  each  collective 
angle  from  5°  to  30°,  the  trend  of  increasing  hub  rotation  rate  could  be  described  with 
a  quadratic  polynomial  fit.  The  quadratic  trends  for  each  collective  angle  fit  the  data 
with  a  least  squares  fit  between  0.98  and  0.999,  except  for  0°  and  -10°  collective. 
The  quadratic  trend  does  not  fit  the  highest  flapping  frequencies  for  0°  collective  and 
-10°  collective,  because  the  hub  rotation  slowed  with  increasing  flapping  frequency  at 
those  points.  For  all  other  data  points,  it  was  a  good  fit. 


♦  0  ■  5  a  10  x  15  x  20  •  25  +30  --10 


Figure  50:  Hub  rotation  vs.  flapping  frequency  for  different  collectives 
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The  absolute  thrust  of  the  Flotor  is  shown  in  Figure  5 1  in  terms  of  blade 
loading.  The  thrust  of  the  blades  is  highest  at  the  maximum  flapping  frequency  and 
25°  collective,  even  though  the  rotation  speed  at  25°  collective  is  not  the  fastest.  At 
25°  collective  and  41  Hz  flapping,  the  blades  produced  18%  more  thrust  despite 
rotating  30%  slower  than  at  15°  collective  and  41  Hz  flapping.  During  steady 
rotation  without  flapping,  the  rotor  thrust  varies  with  the  square  of  rotation  speed  and 
is  approximately  proportional  to  blade  collective  angle.  The  rotation  speed  at  41  Hz 
flapping  and  15°  collective  was  630  RPM.  At  41  Hz  and  25°  collective,  the  rotation 
speed  was  443  RPM.  Rotating  steadily  without  flapping  at  these  two  respective 
speeds,  the  thrust  at  25°  collective  should  be  less  than  the  thrust  at  15°  collective  if  a 
linear  lift  curve  slope  is  assumed.  Under  the  conditions  in  this  test,  however,  the 
thrust  at  25°  collective  is  much  greater,  which  means  that  much  of  the  thrust  at  high 
flapping  frequency  and  high  collective  has  its  source  in  blade  flapping. 

♦  Odeg  ■  5  deg  a  10  deg  x  15  deg 

x  20  deg  •  25  deg  +  30  deg 


Figure  51:  Blade  loading  vs.  flapping  frequency  for  different  collectives 
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The  thrust  coefficient  versus  flapping  frequency  for  all  of  the  collective  angles 
is  shown  in  Figure  52.  The  thrust  coefficient  trends  for  5°,  10°,  and  15°  collective 
show  a  decrease  with  increasing  flapping  frequency.  This  points  to  a  decreasing 
fraction  of  the  thrust  coming  from  flapping  as  flapping  frequency  is  increased.  At 
low  collectives  and  higher  flapping  frequencies,  much  of  the  flapping  energy  is  going 
towards  the  rotation  of  the  blades.  At  higher  collective  angles,  however,  the  trend  of 
thrust  coefficient  remains  constant  or  even  increases  with  flapping  frequency.  With 
increasing  flapping  frequency  at  higher  collectives,  an  increasing  fraction  of  the  thrust 
can  be  attributed  to  blade  flapping. 


0  deg  ■  5  deg  —a—  1 0  deg  x  1 5  deg 


20  deg  —*—25  deg  — i — 30  deg - 10  deg 


Figure  52:  Thrust  coefficient  vs.  flapping  frequency  at  different  collectives 


4.2  Conventional  Powered  Rotation  and  Model  Validation 

The  objective  of  the  next  set  of  tests  was  to  evaluate  the  system  efficiency  of 
the  Flotor  operating  as  a  conventional  MAY  rotor  and  compare  the  results  to  other 
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published  micro  rotor  data.  The  first  experiment  was  conducted  utilizing  powered 
shaft  rotation.  The  flapping  motor  was  not  used  for  these  tests.  The  rotation  motor  is 
an  Astro  Flight  Firefly  coreless  DC  motor  with  a  4: 1  planetary  gearhead.  The  output 
shaft  from  the  motor  gearhead  drives  a  1 : 1  gear,  which  drives  the  rotor  shaft.  The 
motor  speed  was  monitored  with  the  same  hall  sensor  described  earlier,  which  is 
mounted  below  the  rotation  motor  output  gear  in  Figure  53.  The  speed  of  the  motor 
is  manually  controlled  by  varying  the  motor  voltage  from  a  Sorensen  power  supply. 
The  power  supply  has  a  data  output  for  current  and  voltage  which  is  fed  to  the  data 
acquisition  system  in  order  to  monitor  the  electrical  power  consumed  by  the  rotation 
motor.  There  was  not  yet  a  mechanism  in  place  for  measuring  mechanical  power,  so 
any  references  to  power  or  figure  of  merit  in  the  results  are  based  on  electrical  power. 


1:1  Gears 


Magnet 


4:1  Gearhead 


Hall  Sensor 


Firefly  Motor 


Figure  53:  Rotation  motor  mounted  with  hall  sensor 
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The  blades  used  for  this  test  were  carbon  plates  with  a  thickness  ratio  of  3%, 
which  were  molded  with  a  camber  of  about  7%.  Their  shape  was  rectangular  with  a 
chord  of  2  cm  and  a  span  of  8.1  cm.  The  rotor  had  a  total  radius  of  10.8  cm.  This 
type  of  blade  is  used  in  flying  rotary  wing  MAVs  because  of  its  good  performance  at 
low  Reynolds  numbers.  The  data  measured  in  this  experiment  was  the  hub  rotation 
speed,  the  electrical  power  input  to  the  system,  and  the  vertical  thrust  output.  The 
variables  were  the  hub  rotation  speed  and  the  collective  angle.  The  collective  angle 
was  varied  in  steps  of  5°  between  0°  and  25°.  The  hub  rotation  was  tested  from  600 
to  3000  RPM  in  steps  of  200  RPM,  unless  the  electrical  power  capabilities  of  the 
motor  were  exceeded.  For  short  periods,  the  motor  was  tested  beyond  its 
manufacturer  recommended  power  limit  of  10  Watts. 

The  results  for  thrust  versus  RPM  are  shown  in  Figure  54.  The  maximum 
recommended  power  of  the  motor  is  shown  as  a  dotted  line.  As  expected,  the  value 
of  thrust  increases  approximately  with  the  square  of  the  rotor  RPM.  There  is  a  small 
positive  thrust  at  0°  collective  because  the  blades  are  cambered  and  the  airfoil 
sections  produce  lift  even  at  0°  angle  of  attack.  The  greatest  thrust  occurred  at  15° 
collective,  but  exceeded  the  power  recommendations  of  the  motor.  If  constrained  to 
10  Watts  of  power  into  the  motor,  the  thrust  was  approximately  55  grams  at  10°,  15°, 
and  20°  collective  while  rotating  at  2900  RPM,  2400  RPM,  and  2200  RPM, 
respectively. 

Figure  55  shows  the  measured  electrical  power  compared  to  the  rotor  thrust  at 
all  of  the  collective  angles.  The  collective  angles  of  10°,  15°,  and  20°  all  follow  a 
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very  similar  trend  and  all  have  a  similar  maximum  thrust  that  is  limited  by  the  power 
constraint  of  the  motor. 


0  500  1000  1500  2000  2500  3000  3500 

RPM 

Figure  54:  Thrust  vs.  RPM  and  approximate  10  Watt  electrical  power  limit 

At  high  thrust  levels,  the  10°  collective  case  had  the  highest  electrical  power 
loading,  since  it  required  the  least  power  for  a  given  amount  of  thrust.  To  produce  54 
grams  of  thrust  at  10°  collective,  the  motor  required  8.8  Watts  of  power.  To  produce 
the  same  thrust  at  15°  and  20°  collective,  however,  the  motor  required  9.6  W  and  10.3 
W,  respectively. 

The  trend  of  thrust  coefficient  versus  RPM  is  shown  in  Figure  56.  There  is  a 
very  slight  increase  in  thrust  coefficient  with  increasing  RPM.  This  can  be  attributed 
mainly  to  Reynolds  number  effects,  such  as  an  increase  in  lift  curve  slope  and  a 
decrease  in  viscous  effects  with  increasing  rotor  speed.  For  the  most  part,  the  thrust 
coefficient  does  not  vary  appreciably  with  RPM,  therefore,  the  rotor  speed  should 
also  not  have  an  obvious  effect  on  the  aerodynamic  figure  of  merit. 
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Figure  55:  Electrical  Power  vs.  Thrust  for  various  collectives 
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Figure  56:  Ct  vs.  RPM  for  collectives  between  0°  and  25° 

The  figure  of  merit  measured  in  this  experiment  is  the  total  system  figure  of 
merit  based  on  electrical  power  input,  which  is  the  ratio  of  ideal  power  to  actual 
electrical  power.  The  electrical  power  input  into  the  motor  includes  all  losses  such  as 
motor  efficiency,  transmission  efficiency,  and  aerodynamic  efficiency.  The  ideal 
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power  is  the  ideal  aerodynamic  power  based  on  the  measured  rotor  thrust  coefficient. 
At  this  scale,  electric  motors  can  have  a  maximum  efficiency  of  about  70%,  but  only 
when  operated  at  an  ideal  torque  and  RPM.  Their  efficiency  under  normal  operating 
conditions  is  usually  less.  Their  transmissions  are  about  80%  efficient  at  this  size, 
depending  on  design  and  reduction  ratio  [62],  Aerodynamically,  a  figure  of  merit  of 
0.6  is  considered  good  for  an  MAV  rotor  [22,  26],  In  this  case,  the  ratio  of  electrical 
power  to  ideal  power  does  not  exceed  0.2.  Figure  of  merit  is  plotted  versus  RPM  in 
Figure  57. 


Odeg  — ■— 5deg  — lOdeg  x  15deg  — 20deg  — •—  25deg 


For  a  given  collective,  it  has  already  been  shown  that  the  thrust  coefficient 
remains  relatively  constant  with  changes  in  RPM.  As  a  result,  the  trends  for  each 
collective  angle  in  Figure  57  can  be  assumed  to  be  trends  of  constant  thrust 
coefficient.  The  change  in  figure  of  merit  with  RPM  is  therefore  not  due  to  a  change 
in  thrust  coefficient.  The  change  in  figure  of  merit  with  RPM  can  only  be  attributed 
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to  the  changing  mechanical  and  electrical  efficiency  of  the  entire  system.  This 
particular  MAV  system  is  most  efficient  when  rotating  between  1200  and  1500  RPM 
with  the  blades  at  a  collective  angle  of  10°,  which  produces  a  thrust  of  between  9  and 
13  grams.  In  this  range,  the  figure  of  merit  reaches  a  maximum  value  of  0.195.  At 
15°  collective,  the  figure  of  merit  reaches  a  maximum  of  0.19  when  the  blades  are 
rotating  between  1000  and  1200  RPM,  which  also  produces  a  thrust  between  9  and  13 
grams.  When  utilized  in  an  actual  MAV,  the  blades  are  usually  rotated  faster  to 
produce  as  much  thrust  as  possible.  At  maximum  power  and  maximum  thrust,  the 
figure  of  merit  falls  to  about  0.14  at  10°  collective,  0.13  at  15°  collective,  and  0.12  at 
20°  collective.  With  10  Watts  of  input  power,  all  three  of  these  collectives  produce  a 
very  similar  thrust  of  about  55  grams  at  2850,  2400,  and  2180  RPM,  respectively. 
The  lower  tip  speeds  may  be  desirable  to  reduce  aerodynamic  noise,  but  require 
greater  torque. 

The  next  set  of  tests  was  done  to  verify  and  compare  conventional  rotation 
results  of  the  Flotor  with  previously  published  MAV  rotor  results.  Since 
aerodynamic  power  measurement  was  part  of  the  previously  published  results,  a 
sensor  was  added  to  the  Flotor  test  stand  to  measure  aerodynamic  torque.  The  sensor 
was  a  Transducer  Techniques  RTS-25  torque  cell  that  was  mounted  directly  to  the 
base  of  the  test  stand  along  the  axis  of  the  rotor  shaft.  The  full  bridge  circuitry  of  the 
torque  cell  is  easily  integrated  into  the  data  acquisition  system  through  a  SCC-SG04 
module,  which  incorporates  its  own  excitation  voltage  and  signal  conditioning 
amplification. 

The  measurements  for  this  test  included  hub  rotation  speed,  vertical  thrust, 
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rotor  torque,  and  voltage  and  current  into  the  rotation  motor.  The  variables  for  this 
test  were  RPM  and  blade  collective  angle.  The  blade  collective  was  varied  from  2°  to 
36°  in  steps  of  4°.  The  collective  sweep  was  completed  at  660,  1000,  and  1500  RPM. 
The  thrust  data  was  again  gathered  with  the  digital  balance.  The  other  measurements 
were  collected  with  the  data  acquisition  system,  which  collected  data  at  a  rate  of  1000 
points  per  second  for  two  seconds  and  stored  the  average  values  to  an  output  file. 
Each  value  shown  in  a  plot  is  an  average  of  at  least  three  of  the  2-second  averaged 
data  points. 

The  blades  used  for  this  trial  were  previously  tested  by  Hein  and  Chopra  [26] 
to  perform  well  in  low  Reynolds  number  rotor  tests.  The  blade  cross  section  can  be 
seen  in  Figure  58.  They  were  3%  thick  aluminum  plates  with  7%  camber.  Their 
leading  edge  was  sharpened  at  an  angle  of  15°  down  with  respect  to  the  upper  surface. 
The  blades  have  a  rectangular  planform  with  a  chord  of  2  cm  and  a  span  of  6.2  cm. 
The  total  rotor  radius  is  8.95  cm,  which  sweeps  a  total  rotor  area  of  25.2  cm2. 


Leading 

Edge 


3% 


2cm 


Figure  58:  Cambered  plate  airfoil  with  sharpened  leading  edge 


The  results  from  the  experiment  are  summarized  in  Figure  59.  The  previous 
test  is  from  Hein  and  Chopra,  who  tested  the  same  blades  on  a  slightly  smaller  rotor 
with  a  15.2  cm  diameter.  The  tip  Reynolds  number  of  the  Hein  and  Chopra  test  was 
on  the  order  of  40,000,  while  the  maximum  tip  Reynolds  number  during  this  test  was 
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about  20,000  at  1500RPM.  The  thrust  coefficients  during  these  tests  were  at  least 
twice  as  high  as  full  scale  helicopters.  Part  of  the  reason  for  the  high  thrust 
coefficients  is  that  the  solidity  of  these  micro  rotors  is  twice  as  high  as  many  full  scale 
helicopters.  The  effect  of  solidity  is  negated  by  comparing  these  rotors  at  different 
blade  loading  coefficients,  or  CT/o. 


-0—1500  rpm  mech  — © — 1000  rpm  mech  x  660  rpm  mech 

-H — 1500  rpm  elec  — • — 1000  rpm  elec  a  660  rpm  elec 

—  Hein  &  Chopra 


Figure  59:  Mechanical  and  electrical  figure  of  merit  vs  Cj/a 
compared  to  previous  experiment 

At  blade  loading  coefficients  above  0.05,  the  difference  in  performance  with 
varying  RPM  is  obvious.  The  maximum  figure  of  merit  at  660  RPM  is  0.44.  The 
maximum  figure  of  merit  at  1000  RPM  is  0.47  and  the  maximum  in  these  tests  is  0.48 
at  1500  RPM.  The  maximum  figure  of  merit  reached  in  the  tests  by  Hein  and  Chopra 
was  0.52.  The  lower  figure  of  merit  in  these  tests  can  be  attributed  to  a  couple  of 
causes.  First,  the  Reynolds  number  range  of  the  current  tests  is  only  %  to  54  of  the 
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Reynolds  number  of  the  tests  from  Hein  and  Chopra.  This  decrease  in  Reynolds 
number  causes  lower  maximum  lift  coefficients  and  higher  drag  on  the  blades. 
Secondly,  the  rotor  hub  and  blade  grips  in  these  tests  were  designed  as  part  of  a 
robust  flapping  mechanism  and  they  were  not  as  streamlined  as  the  rotor  hub  utilized 
by  Hein  and  Chopra.  Although  all  instruments  were  carefully  calibrated  and  are 
accurate  to  within  1%,  there  are  still  differences  in  calibration  constants  that  could 
account  for  a  small  difference  between  the  two  tests.  The  electrical  figure  of  merit  in 
these  tests  showed  a  slightly  different  trend  than  the  aerodynamic  figure  of  merit. 
The  maximum  electrical  figure  of  merit  was  shifted  to  slightly  higher  blade  loading 
coefficients.  The  maximum  electrical  figure  of  merit  of  0.23  occurred  at  1500  RPM 
over  a  range  of  blade  loading  coefficients  from  0.13  to  0.16.  At  1000  RPM,  the 
electrical  figure  of  merit  peaks  at  a  value  of  0.21  at  a  slightly  higher  blade  loading 
coefficient  range  from  0.16  to  0.19.  The  electrical  figure  of  merit  at  660  RPM  has  the 
lowest  peak  value  of  0.18,  but  peaks  at  the  highest  range  of  blade  loadings  from  0.19 
to  0.21. 

The  experimental  rotor  data  can  now  be  compared  to  the  MATLAB  code 
described  earlier  for  predicting  micro  rotor  performance.  The  blade  element 
momentum  code  is  based  on  lift  and  drag  data  for  a  thin  cambered  plate  without  a 
sharpened  leading  edge.  This  airfoil  data  and  piecewise  fit  used  by  the  code  was 
presented  earlier  in  Figure  33.  The  airfoil  for  the  blades  in  the  current  experiment 
was  similar,  but  the  current  blades  have  a  sharpened  leading  edge.  The  sharpened 
leading  edge  may  cause  slightly  less  drag  than  the  code  predicts. 

The  experimental  and  predicted  thrust  versus  collective  angle  is  shown  in 
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Figure  60.  The  rotor  spinning  at  660  RPM  produces  a  slightly  lower  thrust 
coefficient  at  all  collectives  than  when  spinning  at  1000  or  1500  RPM.  The  code 
prediction  fits  the  experimental  data  extremely  well  at  660  RPM  rotation  and  below 
20  degrees  collective,  possibly  because  the  660  RPM  test  is  closest  in  Reynolds 
number  to  the  airfoil  data  used  in  the  prediction. 


In  the  experiment,  the  thrust  coefficient  increases  linearly  with  increasing 
collective  angle  up  to  22°.  The  model  predicts  the  thrust  coefficient  will  no  longer 
increase  at  a  constant  rate  beyond  21°  collective,  less  than  a  5%  difference  from  the 
experiment.  The  model,  however,  predicts  a  maximum  thrust  coefficient  that  is  15% 
less  than  the  experimental  data  at  a  similar  Reynolds  number.  The  Flotor  operating  as 
a  conventional  rotor  achieved  higher  thrust  coefficients  at  higher  collectives  than 
predicted  by  the  BEMT  model.  Similar  findings  were  discovered  by  Bohorquez  [4] 
with  a  combined  CFD  and  BEMT  model.  He  showed  experimentally  that  a  MAY 
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rotor  stalled  more  gradually  and  reached  thrust  coefficients  up  to  10%  higher  than 
predicted  by  a  BEMT  model  utilizing  CFD  generated  airfoil  data. 

Figure  61  shows  the  experimental  data  and  the  predicted  trend  for  figure  of 
merit  versus  blade  loading  coefficient.  Compared  to  the  data  at  1000  and  1500  RPM, 
the  model  underpredicted  the  figure  of  merit.  Compared  to  the  data  at  660  RPM, 
however,  the  model  predicted  the  figure  of  merit  well  up  to  a  blade  loading 
coefficient  of  0.17.  The  660  RPM  experiment  showed  a  maximum  figure  of  merit  of 
0.44  at  0.19  blade  loading  coefficient.  The  model  predicted  a  maximum  figure  of 
merit  of  0.43  at  a  blade  loading  coefficient  of  0.17,  which  was  only  a  2%  difference 
from  the  experiment.  Above  0.17  blade  loading  coefficient,  the  model  predicted  a 
steep  drop  in  figure  of  merit  that  did  not  occur  in  the  experiment.  These  differences 
are  attributed  to  factors  that  include  the  lower  Reynolds  numbers  of  the  model  data, 
the  sharpened  leading  edge  of  the  experimental  blade,  and  the  inability  for  the  2-D 
BEMT  model  to  accurately  predict  rotor  performance  at  high  collective  angles. 
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Three  versions  of  the  prediction  for  blade  loading  coefficient  versus  collective 
are  shown  in  Figure  62  with  the  660  RPM  experimental  data.  The  first  prediction  is 
the  complete  BEMT  model,  including  a  tip  loss  factor  and  swirling  inflow.  The 
second  prediction  is  the  BEMT  model  including  tip  losses,  but  neglecting  the  swirl 
velocity  component  of  the  inflow.  The  third  prediction  is  the  BEMT  model  without 
tip  losses  or  a  swirl  velocity.  The  best  fit  for  the  data  below  22°  collective  angle  is 
the  complete  model,  which  incorporates  the  swirl  component  of  the  inflow  and  the 
Prandtl  model  for  tip  losses.  Ignoring  the  swirl  component  of  the  inflow  causes  a 
1 0%  increase  in  the  slope  of  the  trend.  Ignoring  both  the  swirl  velocity  and  tip  losses 
in  the  model  causes  a  20  %  increase  in  slope,  a  higher  peak  thrust  prediction,  and  a 
more  abrupt  stall  at  a  lower  collective  angle. 
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Figure  62:  Cx/cr  vs.  collective  for  experimental  data 
and  prediction  with  and  without  losses 
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4.3  Powered  Rotation  with  Powered  Flapping 


The  next  tests  of  the  Flotor  involved  powered  rotation  combined  with 
powered  blade  flapping  at  high  frequencies.  During  the  initial  testing,  measurements 
were  taken  of  average  thrust,  flapping  frequency,  hub  RPM  and  rotation  motor 
electrical  power.  Motor  electrical  power  data  was  input  to  the  data  acquisition  system 
from  a  serial  connection  on  the  back  of  the  Sorensen  power  supply.  The  rotational 
speeds  were  measured  with  Hall  effect  sensors  monitoring  magnets  mounted  on  each 
of  the  two  motor  shafts.  Thrust  data  was  gathered  with  an  Ohaus  Explorer  digital 
balance.  All  of  these  data,  except  for  the  thrust,  were  collected  with  the  National 
Instruments  hardware  and  software  at  a  sample  rate  of  1000  points/second  and  then 
averaged  over  a  period  of  two  seconds.  The  thrust  was  fed  through  a  low  pass  filter 
and  averaged  over  a  period  of  five  seconds  using  the  force  averaging  feature  of  the 
digital  balance.  Mechanical  torque  data  was  not  yet  available  during  these  tests. 

The  first  blades  tested  were  flat  rectangular  aluminum  plates  with  a  thickness 
ratio  of  3%.  They  had  a  chord  of  2  cm  and  a  span  of  6.2  cm.  The  total  rotor  radius 
was  8.95  cm,  which  sweeps  a  total  rotor  area  of  25.2  cm2  with  a  solidity  of  0.14. 
The  blades  were  rotated  at  660  RPM  and  a  tip  speed  of  6.2  m/s.  The  blade  collective 
was  varied  from  0°  to  40°  in  steps  of  4°.  The  low  rotational  speed  was  chosen  for  two 
reasons.  First,  it  closely  represented  the  airspeed  of  MAY  sized  flying  creatures. 
Second,  it  allowed  higher  reduced  frequencies  with  the  same  absolute  frequency  of 
blade  flapping.  The  blade  flapping  frequency  was  4  per  rotor  revolution,  which 
corresponded  to  an  absolute  frequency  of  45  Hz.  The  reduced  frequency  for  this  test 
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was  approximately  0.3  at  the  %  span  of  the  blades.  The  amplitude  of  blade  flapping 
was  ±3°,  which  causes  a  maximum  heaving  velocity  at  the  blade  tips  of 
approximately  2.6  m/s.  Because  of  blade  deformation  during  flapping,  the  angle  of 
attack  variation  could  not  be  accurately  predicted. 

For  each  RPM  and  collective,  the  blades  were  first  rotated  conventionally  at 
the  desired  RPM  and  then  flapped  at  the  desired  frequency.  One  data  point  was 
stored  for  each  case.  A  third  reading  was  recorded  during  zero  rotor  motion  as  a  tare 
test.  Any  residual  values  from  the  tare  test  were  subtracted  from  the  values  gathered 
while  the  rotor  was  in  motion.  Each  data  point  presented  in  the  following  figure 
represents  an  average  of  four  data  points  gathered  under  the  same  conditions. 

The  test  results  are  shown  in  Figure  63.  With  steady  rotation,  the  blade 
loading  coefficient  reached  a  maximum  of  0.16  at  36°  collective,  but  the  thrust  did 
not  increase  much  beyond  28°  collective,  where  the  value  was  0.15.  At  every 
collective  angle,  the  blades  produced  greater  thrust  with  4/rev  flapping  than  they  did 
with  steady  rotation.  Like  the  results  in  the  free  rotation  tests,  the  flexible  blades  in 
this  test  are  producing  extra  thrust  due  to  flapping  at  all  collectives.  The  additional 
thrust  at  collectives  beyond  24°  shows  evidence  of  delayed  stall.  The  thrust  curve 
during  steady  rotation  levels  off  at  collectives  above  24°.  The  thrust  curve  for  the  4 
per  revolution  flapping  case,  however,  maintains  an  almost  linear  slope  up  to  34° 
collective.  This  produces  a  maximum  blade  loading  coefficient  of  0.20  at  36° 
collective,  a  20%  maximum  thrust  increase  over  the  steady  rotation  case. 


94 


0.25 


0.2  • 

•  —  • 

0.15  •  

c  x 

O  •  * 

0.1  #  x 

•  x  x  no  flapping 

0.05  A!  fl  • 

x  •  4/rev  flapping 

0  -f - * - T - T - T - T - 

0  10  20  30  40  50 

collective  (degrees) 

Figure  63:  Blade  loading  coefficient  vs.  collective 
with  660  RPM  rotation  and  4/rev  flapping 

In  order  to  isolate  flexible  blade  phenomena  such  as  rotational  lift  or  purely 
rigid  blade  mechanisms  such  as  dynamic  stall,  the  remainder  of  the  testing  utilized  a 
more  rigid  blade  design.  The  blades  were  the  same  cambered  aluminum  plates  with 
sharpened  leading  edge  that  were  tested  in  the  conventional  rotation  experiment. 
They  have  a  much  higher  area  moment  of  inertia  than  a  flat  plate,  making  them  more 
resistant  to  bending  and  torsion  than  the  flat  plate  blades  in  the  previous  test.  In  the 
remainder  of  the  tests,  the  blades  were  rotated  at  660  and  1000  RPM  and  the 
collective  angle  was  varied  from  0°  to  36°  in  steps  of  2°  or  4°  degrees.  The  blades 
were  tested  with  steady  rotation  and  with  rotation  plus  a  flapping  frequency  between 
2  and  8  flaps  per  rotor  revolution.  The  matrix  of  test  conditions  is  shown  in  Table  1. 
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Figure  63:  Blade  loading  coefficient  vs.  collective 
with  660  RPM  rotation  and  4/rev  flapping 
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Flap  frequency  and  (amplitude) 


None 

2/rev 

4/rev 

6/rev 

8/rev 

(0°) 

(3°) 

(1.5°) 

(1°) 

(0.75°) 

Rotor 

660  RPM 

X 

X 

X 

X 

X 

Speed 

1000  RPM 

X 

X 

X 

X 

Table  1:  Test  matrix  for  rotating  and  flapping  cambered  plate  blades 

Since  the  blades  in  the  following  tests  were  relatively  stiff,  their  angle  of 
incidence  variation  could  be  estimated  with  rigid  blade  element  kinematic  models 
described  in  Chapter  3.  The  blade  angle  of  incidence  is  approximately 

^max  =  ±  tan  1  —  SQ^max  •  A  goal  was  to  keep  the  angle  of  incidence  variation  at  a 

constant  value,  whether  the  blades  were  flapping  at  2/rev  or  8/rev.  Limiting  the 
maximum  angle  of  incidence  during  flapping  keeps  the  blades  from  grossly 
exceeding  stall  angle  of  attack  values  while  at  low  collectives.  The  only  time  the 
blades  should  flap  beyond  their  stall  angle  of  attack  is  when  they  are  set  at  higher 
collective  angles  of  above  20°. 

The  maximum  angle  of  incidence  depends  on  flapping  amplitude,  flapping 
frequency,  and  rotation  speed.  Since  the  rotation  speed  and  flapping  frequencies  for 
these  tests  are  already  constrained,  the  desired  flapping  angle  can  be  calculated. 
Since  the  flapping  will  be  small  angles  in  these  tests,  it  can  be  assumed  that: 
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sinAnax  =  Anax 
tan”*  fcl/7^max  ~  ^fP max 

Q  _  Q 
therefore,  /?max 

Figure  64  shows  the  relationship  between  flapping  angle  and  flapping 
frequency  for  three  different  maximum  angles  of  incidence.  The  goal  for  angle  of 
incidence  variation  for  these  tests  was  chosen  as  approximately  ±  6°.  The  four  points 
on  the  line  for  6°  angle  of  incidence  correspond  to  the  four  flapping  frequencies 
presented  in  the  next  set  of  tests.  The  flapping  amplitude,  /3max  is  varied  by  altering 
the  crank  radius  of  the  flywheel  on  the  flapping  motor.  The  crank  radius  remained 
the  same  for  the  tests  of  660  and  1000  RPM  at  the  same  number  of  flaps  per  rotor 
revolution.  The  crank  radius  was  altered  for  each  test  of  a  different  number  of  flaps 
per  revolution.  It  will  be  shown  later  that  the  actual  flapping  amplitude  may  have 
differed  from  the  calculated  values  due  to  factors  such  as  part  tolerances  and  inertial 
effects  of  the  flapping  mechanism  at  high  frequencies.  Dynamic  measurement  of  the 
flapping  amplitude  will  be  described  with  the  results. 

In  the  remainder  of  Flotor  tests,  measurements  were  taken  of  average  thrust, 
average  rotor  torque,  flapping  motor  RPM,  rotation  motor  RPM  and  rotation  motor 
electrical  power.  All  data  were  again  gathered  and  averaged  with  the  data  acquisition 
system  except  for  the  thrust,  which  was  averaged  with  the  digital  balance.  Data  were 
gathered  during  each  of  the  three  conditions  of  conventional  rotation,  rotation  plus 
flapping,  and  zero  motion.  Each  data  point  shown  for  each  test  condition  is  an 
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average  of  at  least  three  data  samples  gathered  under  the  same  conditions. 


Figure  64:  Flapping  amplitude  vs.  flapping  frequency  ratio 
for  constant  values  of  angle  of  incidence 


The  data  for  each  RPM  is  presented  in  terms  of  blade  loading  coefficient 
versus  collective,  torque  coefficient  versus  blade  loading  coefficient,  and  figure  of 
merit  versus  blade  loading  coefficient.  The  figure  of  merit  is  more  for  informational 
purposes,  since  it  does  not  include  the  power  input  to  the  flapping  motor.  Since  the 
blade  loading  is  directly  proportional  to  the  thrust,  sometimes  a  change  in  blade 
loading  coefficient  will  be  described  as  a  change  in  thrust.  Similarly,  a  change  in 
torque  coefficient  is  sometimes  referred  to  as  a  change  in  torque. 

In  the  three  samples  of  the  same  test  condition,  the  standard  deviation  is 
typically  1%  for  thrust,  3%  for  torque  and  2%  for  RPM.  These  errors  will  propagate 
to  an  expected  error  of  approximately  4%  in  thrust  coefficient  and  7%  in  power 
coefficient,  both  of  which  are  used  in  the  calculation  of  figure  of  merit.  The  figure  of 
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merit  will  have  an  expected  error  of  approximately  8%,  which  is  the  reason  the  figure 
of  merit  data  had  the  most  scatter.  The  following  data  was  previously  presented  by 
Fitchett  and  Chopra  [63],  with  the  exception  of  the  tests  at  1000  RPM  with  6  per 
revolution  flapping  and  660  RPM  with  8  per  revolution  flapping. 

The  first  test  incorporated  2/rev  flapping  frequency  at  an  amplitude  of  3°.  The 
blades  were  rotated  at  660  and  1000  RPM  over  a  range  of  blade  collectives.  The 
reduced  frequency  at  the  3A  blade  span  during  flapping  and  rotation  was 
approximately  0.15.  The  data  with  flapping  and  without  flapping  has  been  plotted 
simultaneously  with  different  symbols.  The  660  RPM  data  for  this  test  is  shown  in 
Figure  65  and  the  1000  RPM  data  is  shown  in  Figure  66.  The  blade  loading 
coefficient  for  660  RPM  reaches  a  maximum  of  0.21  at  30°  collective  with  no 
difference  between  the  flapping  and  non-flapping  case.  The  blade  loading  coefficient 
in  the  1000  RPM  case  reaches  a  maximum  of  0.22  at  32°  collective  for  both  the 
flapping  and  non  flapping  case.  At  high  collective  angles  of  34°  and  36°,  the  1000 
RPM  2/rev  flapping  case  shows  an  increase  in  thrust  of  about  5%  compared  to  the 
non-flapping  case,  but  this  does  not  reach  beyond  the  maximum  thrust  of  either  case. 

The  torque  coefficients  for  660  and  1000  RPM  with  steady  rotation  are  within 
5%  of  each  other  at  every  collective  angle.  The  only  differences  occur  with  flapping 
at  high  collectives.  Even  though  the  1000  RPM  case  achieves  slightly  increased 
thrust  with  flapping  at  the  two  highest  collective  angles,  the  torque  at  those  points  is 
almost  identical  to  the  torque  without  flapping.  For  flapping  and  steady  rotation,  the 
figure  of  merit  reaches  a  maximum  of  0.46  at  1000  RPM  and  0.43  at  660  RPM.  In 
general,  there  is  no  appreciable  difference  or  benefit  due  to  2/rev  blade  flapping. 
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Figure  65:  660  RPM  2/rev  flap  results 
(a)  C t/c  vs.  collective,  (b)  Cq  vs.  Cj/g  (c)  FM  vs.  Cj/g 
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Figure  66:  1000  RPM  2/rev  flap  results 
(a)  CT/o  vs.  collective,  (b)  Cq  vs.  Cj/g  (c)  FM  vs.  C  |/ct 
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The  next  test  was  conducted  with  the  same  rotational  speeds,  but  the  flapping 
frequency  was  increased  to  4/rev  at  an  amplitude  of  1.5°.  During  blade  flapping,  the 
reduced  frequency  was  0.3  at  %  rotor  radius.  The  660  RPM  data  are  shown  in  Figure 
67  and  the  1000  RPM  data  are  shown  in  Figure  68.  Above  26°  collective  and  below 
8°  collective,  there  was  an  increase  in  thrust  with  4/rev  flapping.  With  flapping,  the 
maximum  thrust  increased  by  6%  at  660  RPM  and  32°  collective  and  similarly  by  5% 
at  1000  RPM  and  30°  collective.  Beginning  to  show  here  is  evidence  of  delayed  stall 
at  collective  angles  above  24°.  During  the  cases  with  4/rev  flapping,  the  blades  were 
actually  stalling  at  slightly  higher  collective  angles  than  in  the  cases  without  flapping. 
The  linear  slope  of  the  thrust  versus  collective  plot  extends  to  collective  angles  about 
4°  higher. 

The  increased  maximum  figure  of  merit  for  the  4/rev  flapping  cases  is  not 
only  attributed  to  the  increased  thrust,  but  also  to  the  reduction  of  torque  necessary  to 
spin  the  rotor.  With  flapping,  the  torque  was  reduced  by  10%  when  operating  at 
maximum  FM,  which  occurred  at  a  Cj/o  between  0.1  and  0.15.  At  moderate 
collectives,  the  torque  reduction  was  evidence  that  the  blades  were  producing  some  of 
their  own  force  in  the  direction  of  rotation.  At  collective  angles  above  30°,  the  torque 
during  flapping  increased  by  4  -  6%,  showing  signs  of  increased  drag  on  the  blades. 
Okamoto  and  Azuma  [37]  also  measured  similar  drag  increases  during  dynamic  stall 
of  heaving  wings  at  similar  Reynolds  number  and  reduced  frequency. 
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Figure  67:  660  RPM  4/rev  flap  results 
(a)  CT/o  vs.  collective,  (b)  Cq  vs.  Cj/g  (c)  FM  vs.  C  |/ct 
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Figure  68:  1000  RPM  4/rev  flap  results 
(a)  Ct/cj  vs.  collective,  (b)  Cq  vs.  C j/a  (c)  FM  vs.  Cj/a 
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The  next  test  was  conducted  at  660  RPM  with  6/rev  flapping  at  a  target 
flapping  amplitude  of  1°.  It  should  be  noted  that  the  actual  flapping  amplitude  was 
nearly  2°  during  this  test  due  to  tolerances  in  the  flapping  mechanism.  Reasons  for 
the  discrepancy  and  a  method  for  measuring  flapping  angle  will  be  explained  in  detail 
after  all  of  the  results  are  presented.  Nonetheless,  the  6/rev  reduced  frequency  met 
the  goal  of  0.45  since  it  does  not  depend  on  flapping  amplitude.  The  data  for  this  test 
are  in  Figure  69. 

The  thrust  coefficient  at  all  collectives  increased  with  6/rev  flapping  compared 
with  conventional  rotation.  In  this  case,  it  is  apparent  that  delayed  stall  was 
occurring,  which  allowed  the  rotor  with  6/rev  flapping  to  achieve  average  maximum 
thrust  coefficients  15%  higher  than  conventional  rotation  at  30°  collective.  The 
flapping  blades  were  again  producing  a  force  in  the  direction  of  rotation,  evidenced 
by  an  obvious  reduction  in  the  shaft  torque  below  30°  collective.  With  6/rev 
flapping,  the  torque  was  reduced  by  at  least  30%  at  collective  angles  less  than  16°  and 
the  torque  was  reduced  by  10%  while  the  rotor  operated  at  maximum  figure  of  merit. 
Torque  was  equal  to  the  conventional  rotation  case  when  operating  at  maximum 
thrust.  The  torque  data  with  6/rev  flapping  showed  a  3%  increase  at  collective  angles 
above  30°. 

The  next  results  in  Figure  70  are  from  the  test  with  1000  RPM  rotation  and 
6/rev  blade  flapping  at  a  target  amplitude  of  1°.  The  reduced  frequency  was  0.45  at  3A 
rotor  radius.  The  error  of  flapping  amplitude  in  this  case  was  within  30%,  which  is 
much  improved  compared  to  previous  test.  The  data  again  show  an  increase  in  thrust 
with  flapping  at  collectives  below  8°.  At  0°  collective,  the  blade  loading  coefficient 
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with  flapping  is  0.026  compared  to  0.007  with  conventional  rotation. 
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Figure  69:  660  RPM  6/rev  flap  results 
(a)  Ct/c  vs.  collective,  (b)  Cq  vs.  Ct/cj  (c)  FM  vs.  Ct/ct 
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The  thrust  data  again  also  shows  evidence  of  dynamic  stall  occurring  at 
collective  angles  above  24°.  With  flapping,  the  linear  slope  of  the  blade  loading 
coefficient  data  extends  4-6°  higher  than  conventional  rotation  before  the  thrust  data 
show  signs  of  blade  stall.  The  maximum  blade  loading  coefficient  with  steady 
rotation  is  0.216  at  30°  collective.  With  flapping,  the  maximum  blade  loading 
coefficient  increases  by  13%  to  0.243  at  34°  and  36°  collective. 

The  torque  data  shows  obvious  effects  of  delayed  stall  at  high  collective 
angles.  With  6/rev  flapping  above  28°  collective,  the  torque  shows  an  8%  to  13% 
increase  compared  to  conventional  rotation.  The  torque  decrease  at  lower  collective 
angles  in  this  test  was  not  as  significant  as  was  seen  in  the  previous  test  of  higher 
amplitude  6/rev  flapping.  The  decrease  in  torque  in  the  previous  tests  at  moderate 
collectives  seems  more  dependent  on  flapping  amplitude  than  it  is  on  flapping 
frequency.  Some  of  this  phenomenon  can  be  explained  quasi-steadily.  On  a  blade 
element  at  radius  y,  the  lift  vector  dL  tilts  forward  by  the  angle  of  incidence  (j>  during 
the  downstroke  and  the  upstroke,  producing  a  blade  element  “anti-torque”  of  ydLs'm<f>. 
If  the  angle  of  incidence  is  increased  by  increasing  the  flapping  amplitude,  then  the 
lift  vector  will  produce  greater  anti  torque.  According  to  the  data,  this  quasi-steady 
explanation  is  either  only  valid  at  higher  absolute  flapping  amplitudes  or  it  starts 
becoming  invalid  at  high  flapping  frequencies. 
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Figure  70:  1000  RPM  6/rev  flap  results 
(a)  Ct/cj  vs.  collective,  (b)  Cq  vs.  Ct/cj  (c)  FM  vs.  CT/a 
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The  final  test  completed  in  this  series  of  rigid  blade  tests  involved  rotation  at 
660  RPM  and  blade  flapping  at  8  per  revolution.  The  reduced  frequency  was  0.6  at  14 
blade  span.  The  flapping  magnitude  target  was  0.75°,  which  was  found  to  be  smaller 
than  the  minimum  possible  flapping  amplitude  of  the  Flotor  at  this  frequency.  The 
minimum  flapping  amplitude  at  high  frequencies  is  apparently  limited  by  the 
tolerances  and  backlash  of  the  components.  The  actual  flapping  amplitude  was 
slightly  greater  than  1°,  which  was  achieved  with  a  flapping  motor  crank  radius  of 
0.007”  ±0.001”.  The  test  results  are  shown  in  Figure  71. 

The  data  below  8°  collective  with  8/rev  flapping  again  shows  an  increase  in 
thrust  compared  to  conventional  rotation,  but  it  is  much  less  obvious  than  the 
previous  4/rev  or  6/rev  tests.  At  0°  collective,  the  blade  loading  coefficient  with 
flapping  is  0.015  compared  to  the  steady  rotation  value  of  0.010.  The  thrust  increase 
at  high  collective  angles  is  more  familiar.  With  8/rev  flapping,  the  maximum  blade 
loading  coefficient  is  0.240  at  34°  collective,  which  is  a  13%  increase  compared  to 
the  conventional  rotation  value  of  0.213  at  30°  collective.  The  torque  at  low  thrust 
values  with  flapping  again  shows  a  slight  decrease  compared  to  conventional  rotation. 
The  torque  at  high  thrust  values  again  shows  evidence  of  dynamic  stall.  With 
flapping,  the  torque  values  at  32°,  34°,  and  36°  collective  increase  by  5%,  10%,  and 
15%,  respectively,  compared  to  conventional  rotation. 
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Figure  71:  660  RPM  8/rev  flap  results 
(a)  Ct/ct  vs.  collective,  (b)  Cq  vs.  Cj/o  (c)  FM  vs.  Cj/a 
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In  the  flapping  results,  especially  at  6/rev,  an  interesting  thrust  increase  occurs 
near  0°  collective.  In  Figure  69(a),  at  0°  collective,  the  rotor  produced  a  blade 
loading  coefficient  with  6/rev  flapping  of  .027  compared  to  .005  without  flapping.  In 
Figure  70(a),  at  0°  collective,  the  rotor  produced  a  blade  loading  coefficient  with 
6/rev  flapping  of  .025  compared  to  .007  without  flapping.  This  thrust  increase  was 
also  apparent  in  the  4/rev  tests,  albeit  not  as  much.  It  can  be  explained  by  purely 
steady  phenomena  due  to  the  variation  of  lift  and  drag  vectors  with  angle  of  attack. 
Recalling  the  lift  and  drag  curves  of  a  cambered  plate  airfoil  from  Figure  33,  the  lift 
curve  slope  is  9  times  greater  above  -4°  angle  of  attack  than  below.  The  drag 
increases  faster  at  positive  angles  of  attack  compared  with  negative  angles  of  attack. 
For  the  purpose  of  this  example,  assume  that  a  blade  section  is  at  0°  collective  and 
undergoes  an  angle  of  attack  variation  of  ±10°  due  to  flapping.  A  detailed  view  of 
this  portion  of  the  lift  and  drag  curves  is  in  Figure  72  along  with  the  maximum  values 
of  lift  and  drag  coefficient  due  to  10°  angle  of  attack  variation. 


Figure  72:  Lift  and  drag  vs.  angle  of  attack  for  cambered  plate  including 
maximum  and  minimum  values  for  ±10°  a  variation 
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When  the  blade  is  flapping  downwards,  the  angle  of  attack  reaches  a 
maximum  of  10°  and  the  drag  coefficient  reaches  a  maximum  of  0.213.  The  drag 
vector  on  the  downstroke  contributes  to  the  vertical  force  of  the  rotor,  as  shown 
previously  in  Figure  34.  The  portion  of  the  drag  coefficient  acting  upwards  during 
downward  flapping  is  0.2 13 sin  10°  =0.037  and  the  vertical  portion  of  drag  during 
upward  flapping  is  0.159sin(-10°)  =  -0.028 .  Since  the  drag  curve  is  not  symmetric 
about  0°  angle  of  attack,  the  drag  averages  to  a  net  upward  force  coefficient  of  0.0045 
while  flapping  at  0°  collective. 

The  vertical  component  of  lift  coefficient  during  downward  flapping  is 
1.13cos(10°)  =  l.ll  and  during  upward  flapping  is  - 0.3 cos(-10°)  =  -0.295.  Also 
because  the  lift  curve  is  not  symmetric  about  0°  angle  of  attack,  the  lift  produces  a  net 
force  upwards  during  blade  flapping.  The  lift  will  produce  an  average  vertical  force 
coefficient  of  0.41  during  flapping,  which  is  more  than  double  the  steady  lift 
coefficient  of  0.17  at  0°  angle  of  attack.  The  difference  in  thrust  seen  while  flapping 
at  low  collective  angles  is  a  quasi-steady  effect  that  is  primarily  due  to  the  properties 
of  the  asymmetric  lift  curve  slope  of  highly  cambered  blades  and  secondarily  due  to 
the  vertical  contribution  of  the  drag. 

As  mentioned  earlier,  the  results  for  660  RPM  6/rev  flapping  were  drastically 
different  from  the  previous  4/rev  results  because  the  blade  flapping  amplitude  goal 
was  being  overshot.  The  backlash  and  tolerances  of  the  parts  were  allowing  the 
inertia  of  the  blades  to  carry  them  to  greater  angles  than  intended.  The  tolerances  of 
the  parts  were  within  0.001”  and  the  bearings  in  the  load  path  are  ABEC-5  rated.  At 
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higher  flapping  frequencies,  however,  the  driving  pin  in  the  flapping  motor  flywheel 
is  offset  by  less  than  0.01”,  which  means  there  is  at  least  10%  error  built  into  the 
flapping  mechanism.  The  backlash  error  percentage  increases  as  the  flapping 
amplitude  decreases. 

The  motion  of  the  test  setup  was  experimentally  investigated  further  and 
corrected  after  the  660  RPM  6/rev  test.  Prior  to  this  point,  the  flapping  amplitude  was 
measured  by  slowly  rotating  the  flapping  motor  and  taking  measurements  of  the  blade 
angle  at  its  peak  values,  which  is  equivalent  to  measuring  the  angles  at  very  low 
flapping  frequencies.  In  order  to  measure  the  flapping  at  higher  target  frequencies, 
the  motion  was  photographically  analyzed.  A  Canon  A530  digital  camera  was 
utilized  in  extended  exposure  mode  to  capture  the  flapping  amplitude  at  high 
frequencies.  The  camera  was  secured  on  a  level  surface  at  the  same  height  as  the 
rotor  plane  and  the  shutter  was  set  to  open  for  at  least  two  seconds  to  capture  the 
maximum  peak  to  peak  blade  flapping  amplitude.  The  photographs  were  then 
analyzed  graphically  with  Engauge  Digitizer,  which  allows  accurate  digital 
measurement  of  photographic  data  or  other  images  [64],  Figure  73  is  a  photo  of  the 
Flotor  flapping  at  about  60  Hz  with  a  target  amplitude  of  2°  peak  to  peak.  The  actual 
measured  peak  to  peak  amplitude  in  this  photo  is  closer  to  5°.  It  was  found  that,  at 
the  2  and  4/rev  flapping  frequencies,  the  flapping  amplitude  was  reaching  its  targets 
within  30%,  but  the  6/rev  flapping  frequency  was  more  than  doubling  the  target 
amplitude.  This  large  uncertainty  with  the  current  mechanism  was  difficult  to  avoid 
due  to  the  tolerances  already  described.  For  the  remainder  of  the  tests,  this  overshoot 
was  corrected  by  reducing  the  crank  radius  of  the  offset  pin  on  the  flywheel 
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accordingly.  By  the  final  test,  the  pin  was  offset  by  only  0.007”. 


Figure  73:  Flotor  blade  with  5°  peak  to  peak  flapping 


To  demonstrate  the  full  effects  of  unsteady  aerodynamics  on  the  average  rotor 
forces,  the  analytical  BEMT  model  is  again  employed.  The  model  calculates  the 
purely  steady  aerodynamic  forces,  so  it  can  be  used  to  observe  which  phenomena  in 
the  experiment  are  due  to  quasi-steady,  two  dimensional  effects.  The  thrust  data  for 
the  quasi-steady  model  are  shown  in  Figure  74. 


Figure  74:  Cj/o  vs.  collective  (quasi  steady  prediction  for  4/rev  flapping) 
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The  quasi-steady  model  predicts  the  increase  in  thrust  coefficient  at  collective 
angles  below  5°.  At  0°  collective,  the  model  predicts  that  the  blade  loading 
coefficient  during  flapping  will  be  0.015  compared  with  0.007  for  steady  rotation. 
Experimental  results  showed  blade  loading  coefficients  of  about  0.25  with  6/rev 
flapping  at  0°  collective.  This  phenomena  is  predicted  by  steady  aerodynamics. 
According  to  the  quasi-steady  model,  the  increased  thrust  at  0°  collective  will 
increase  with  increasing  flapping  amplitude.  The  quasi-steady  prediction  does  not 
predict  well  the  thrust  of  the  flapping  rotor  above  15°  collective  angle.  The  variation 
of  angle  of  attack  is  bringing  the  quasi-steady  model  beyond  stall,  so  the  average  lift 
decreases.  This  is  where  the  unsteady  effects  are  seen  most  in  the  experimental  data. 
If  it  were  not  for  delayed  stall  above  15°  collective,  the  thrust  with  4/rev  flapping 
would  actually  be  decreased  compared  to  conventional  rotation. 

The  quasi  steady  model  of  torque  versus  collective  is  shown  in  Figure  75. 
The  features  of  the  torque  prediction  at  low  to  moderate  collective  angles  match  very 
well  with  the  experimental  data  trends.  At  collectives  below  5°  with  4/rev  flapping, 
there  is  less  torque  required  than  conventional  rotation.  At  collectives  between  5°  and 
15°  with  4/rev  flapping,  there  is  less  torque  required  to  produce  the  same  amount  of 
thrust  as  conventional  rotation.  Above  25°,  the  quasi-steady  model  predicts  the  same 
amount  of  torque  with  less  thrust  for  the  4/rev  flapping  case.  The  experimental  data, 
however,  shows  more  torque  and  more  thrust  at  the  highest  collective  angles  with 
flapping.  The  quasi-steady  prediction  is  adequate,  except  at  high  collectives  with 
high  frequency  flapping.  This  is  the  region  where  the  unsteady  effects  dominate  the 
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forces  on  the  rotor. 


Figure  75:  Torque  vs.  collective  (quasi-steady  prediction  for  4/rev  flapping) 


5  Summary  and  Concluding  Remarks 

The  goal  of  this  study  was  to  introduce  and  examine  the  effects  of  biologically 
inspired  motion  on  current  hover  capable  technology.  After  studying  the  methods 
used  by  nature  for  efficient  flight,  avian  based  kinematics  were  emulated  into  a 
simple  flapping  motion  that  could  be  applied  to  a  current  MAV  rotor.  It  was 
hypothesized  and  proven  that  this  combination  of  rotary  and  flapping  motion  could 
provide  hover  capable  thrust  without  the  necessity  for  powered  shaft  rotation.  This 
thesis  described  the  design,  fabrication  and  testing  of  a  hover  capable  micro  air 
vehicle  rotor  capable  of  being  powered  only  by  avian  based  blade  flapping. 
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5.1  Contributions  and  Conclusions 


The  Flotor  research  test  bed  was  designed  and  manufactured  for  conventional 
and  innovative  flapping  MAV  rotor  experiments.  This  experience  demonstrated  the 
importance  of  careful  part  and  assembly  design.  Design  iterations  were  more  quickly 
performed  on  a  3-D  computer  model  so  that  only  one  complete  hardware  assembly 
was  constructed.  A  digital  3-D  catalog  of  each  part  of  the  assembly  is  available  to  be 
used  by  any  interested  researcher.  The  Flotor  was  capable  of  operating  in  each  of  its 
three  hover  capable  motions:  powered  blade  flapping  with  passive  (unpowered)  hub 
rotation,  conventional  powered  rotation,  and  powered  rotation  with  powered  flapping. 
A  blade  element  momentum  analysis  was  conducted  for  the  cases  of  powered  rotation 
with  and  without  blade  flapping. 

5.1.1  Flapping  blades  with  passive  hub  rotation 

The  first  tests  of  the  Flotor  proved  that  flexible  blades  with  avian  based  blade 
flapping  could  cause  the  blades  to  rotate  with  zero  torque  input  to  the  rotor  shaft. 
These  were  perhaps  the  first  experiments  of  a  torqueless  MAY  flapping  rotor  with  no 
motor  to  power  the  rotor  shaft.  The  first  tests  were  conducted  with  two  sets  of 
flexible  blades,  one  of  which  emulated  the  shape  a  bat’s  wings.  The  results  from 
those  tests  of  flapping  with  passive  rotation  led  to  the  following  conclusions: 

1.  The  hub  rotation  speed  increased  at  a  greater  rate  than  the  blade  flapping 
frequency.  The  hub  rotation  reached  the  fastest  steady  values  at  collective 
angles  of  15°-20°. 

2.  At  positive  collective  angles,  both  the  blade  flapping  and  the  hub  rotation 
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contributed  to  the  vertical  thrust.  The  fraction  of  the  thrust  attributed  to  blade 


flapping  grew  with  increased  collective  angle  and  increased  blade  flapping 
frequency.  The  greatest  thrust  was  a  certain  combination  of  translational  lift 
due  to  hub  rotation  plus  vertical  force  from  flexible  blade  flapping  at  a  high 
collective  angle. 

3.  Maximum  thrusts  and  blade  loadings  of  10  N/m  were  close  to  those  of 
similarly  sized  flying  creatures,  but  they  are  far  below  the  values  for  current 
MAVs. 

5.1.2  Conventional  rotation 

The  second  set  of  experiments  reinforced  the  limited  data  available  for 
conventional  MAY  rotor  performance.  It  also  compared  aerodynamic  efficiency  with 
total  propulsion  system  efficiency  in  terms  of  ideal  aerodynamic  power  and  actual 
electrical  power.  The  Flotor  was  tested  with  conventional  powered  shaft  rotation  and 
two  sets  of  cambered  plate  blades,  one  with  a  sharpened  leading  edge.  The 
conventional  rotation  data  was  used  to  verify  a  BEMT  analysis.  Prandtl  tip  losses  and 
a  swirl  velocity  derived  from  conventional  propeller  theory  were  both  incorporated 
into  the  BEMT  model.  The  following  observations  were  apparent: 

1.  The  operating  conditions  for  maximum  aerodynamic  efficiency  are  not  the 
same  as  the  conditions  for  maximum  propulsion  system  efficiency.  The 
maximum  system  efficiency  occurred  at  higher  collectives  and  thrust  levels 
than  those  for  maximum  aerodynamic  efficiency.  Matching  the  propulsion 
system  to  the  rotor  is  more  important  than  the  maximum  efficiency  of  either 
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component. 

2.  Limited  to  a  maximum  electrical  power,  the  system’s  maximum  thrust 
remained  relatively  constant  between  10°  and  20°  collective  angle.  The 
maximum  thrust  at  20°  collective  required  lower  tip  speeds  and  was 
objectively  quieter  than  at  10°  collective  at  the  same  maximum  thrust. 

3.  The  present  experimental  rotor  data  fits  very  well  with  previously  published 
data  [26],  Because  of  the  lower  tip  speeds  of  the  current  study,  there  were 
some  detrimental  effects  of  reduced  Reynolds  number,  such  as  decreased 
maximum  thrust  coefficients  and  a  figure  of  merit  that  does  not  exceed  0.5.  In 
the  Reynolds  number  range  of  10,000  -  40,000,  doubling  the  Reynolds 
number  made  possible  an  8%  increase  in  maximum  figure  of  merit. 

4.  For  normal  operating  conditions  below  21°  collective,  the  BEMT  analysis 
predicted  the  rotor  performance  and  maximum  figure  of  merit  very  well. 
Above  21°  collective,  however,  the  analysis  predicted  an  earlier  rotor  stall 
than  shown  by  the  experiments.  The  experimental  rotor  achieved  higher 
thrust  at  higher  collective  angles  than  predicted  by  the  BEMT  analysis.  This 
observation  agreed  with  previously  published  comparisons. 

5.1.3  Powered  hub  rotation  with  powered  blade  flapping 

The  final  set  of  results  was  from  obtained  powered  hub  rotation  plus  powered 
blade  flapping.  One  test  was  accomplished  with  flexible  flat  plate  blades  flapping  at 
a  frequency  of  4/rev.  Relatively  rigid  cambered  plates  were  tested  over  a  range  of 
flapping  frequencies  up  to  8  per  rotor  revolution,  or  a  reduced  frequency  of  0.6.  All 
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results  were  simultaneously  compared  to  results  of  conventional  rotation.  The  BEMT 
analysis  developed  earlier  was  used  to  explain  some  features  of  the  results.  The 
results  showed  the  following: 

1.  The  flapping  flexible  flat  plate  blades  produced  a  thrust  benefit  at  every 
collective  angle.  At  high  collectives,  the  flapping  blades  showed  a  large 
benefit,  increasing  the  maximum  thrust  by  20%  due  to  delayed  stall. 

2.  Benefits  of  rigid  cambered  blade  flapping  at  frequencies  greater  than  2/rev 
included  increased  maximum  thrust  coefficients  by  up  to  15%  and  decreased 
torque  required  at  moderate  thrust.  There  was  a  thrust  produced  with  flapping 
near  0°  collective  due  to  the  asymmetric  lift  and  drag  curves  of  the  highly 
cambered  blades.  The  maximum  blade  loading  coefficient  of  0.251  for  the 
cambered  plate  blades  with  6/rev  flapping  was  possibly  the  highest  recorded 
for  a  rotary  wing  MAV. 

3.  At  the  high  flapping  amplitudes  during  6/rev  flapping  at  collective  angles  less 
than  16°,  shaft  torque  was  reduced  by  an  average  of  more  than  30%. 

4.  The  quasi-steady  BEMT  analysis  fit  the  experimental  data  well  below  15° 
collective,  with  and  without  blade  flapping.  The  model  accurately  predicted  a 
reduced  torque  with  blade  flapping  at  moderate  collectives  and  an  increase  in 
thrust  at  0°  collective.  Above  15°  collective,  however,  the  model  could  not 
predict  even  the  average  loads  measured  in  the  experiments. 

5.2  Future  Work  Recommendations 

A  conventionally  powered  rotor  implementing  an  additional  actuation  source 
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for  blade  flapping  is  too  complicated  for  a  flying  vehicle.  In  addition,  blade  flapping 
to  achieve  additional  thrust  from  delayed  stall  requires  high  collective  angles  and 
impairs  rotor  efficiency  due  to  the  increased  shaft  torque  required.  Instead,  flexible 
blade  flapping  at  larger  amplitudes  should  be  implemented  in  order  to  eliminate  the 
torque  required  during  normal  operating  conditions. 

Modeling  of  a  flapping  rotor  with  BEMT  and  steady  aerodynamics  has  been 
adequate  for  initial  modeling  of  the  aerodynamics  of  this  type  of  motion.  A  simple 
elastic  blade  model  should  be  incorporated  to  estimate  the  blade  aerodynamic  angles. 
CFD  and  flow  visualization  may  be  necessary  to  estimate  the  behavior  of  the  time 
varying  inflow. 

To  determine  feasibility  of  a  torqueless  rotor  flying  vehicle,  the  aerodynamic 
performance  must  first  be  improved  in  a  number  of  areas.  The  flexible  blade  design 
used  in  the  torqueless  configuration  needs  to  be  explored.  The  blade  aerodynamic 
properties  should  incorporate  twist,  camber,  anhedral,  and  taper.  They  need  to  be 
aeroelastically  tailored  to  provide  a  balance  of  vertical  force  for  thrust  and  horizontal 
force  to  power  the  rotation  of  the  rotor.  During  the  downstroke,  the  blade  should 
provide  a  large  vertical  force  plus  all  of  the  horizontal  force  necessary  to  power 
rotation  of  the  rotor.  The  beginning  of  the  downstroke  may  consume  the  vast 
majority  of  the  power  of  the  flapping  motion.  During  the  upstroke,  the  blade  should 
harness  the  momentum  of  the  rotor  and  utilize  airspeed  to  produce  lift.  The  lift  will 
assist  in  accelerating  the  blade  upward  at  the  end  of  the  downstroke  and  beginning  of 
the  upstroke.  These  motions  are  modeled  after  a  bird  in  forward  flight.  The  motions 
and  aerodynamics  can  be  examined  with  the  current  Flotor  test  bed. 
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It  has  proven  almost  impossible  with  the  current  instruments  and 
configuration  to  measure  instantaneous  aerodynamic  forces  on  the  blades  of  the 
Flotor  because  reciprocating  forces  of  the  moving  parts  can  be  orders  of  magnitude 
higher  than  the  aerodynamic  forces.  Different  methods  for  finding  the  aerodynamic 
forces  should  be  explored.  The  aerodynamic  forces  could  be  increased  and  the 
reciprocating  mass  forces  could  be  decreased  by  operating  the  flotor  blades  at  a 
slower  flapping  frequency  in  a  more  viscous  fluid.  The  fluid  forces  would  become 
the  dominant  forces  on  the  mechanism  and  the  fluid  flow  could  be  easily  visualized 
and  studied.  Possible  fluids  include  water  or  oil  of  varying  viscosities.  Similar 
methods  have  been  successfully  used  by  many  researchers  studying  insect  flapping 
wing  aerodynamics. 

Once  the  aerodynamics  and  blade  structure  is  optimized,  it  must  then  be 
determined  if  a  flapping  rotor  system  can  be  designed  light  enough  for  a  flying 
vehicle.  Powered  rotary  motion  is  not  necessary  for  blade  flapping,  so  the  idea  of  a 
motor  comprised  of  dense  copper  and  ferrous  materials  should  be  abandoned.  Instead 
of  a  system  that  converts  rotary  to  linear  motion  as  on  the  current  Flotor,  a  blade 
flapping  mechanism  should  be  designed  to  operate  at  resonance.  Alternative  systems 
for  high  frequency  periodic  actuation  should  be  investigated,  such  as  chemical 
muscles  or  smart  materials.  On  this  scale,  utilizing  a  lightweight  and  efficient 
periodic  actuation  source  coupled  with  a  flapping  mechanism  operating  at  resonance 
should  require  drastically  less  power  than  a  conventional  electric  motor.  As  size 
decreases,  this  type  of  system  may  prove  more  efficient  than  the  conventional  scaled 
down  helicopter  rotors  that  are  currently  utilized  for  MAVs. 
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